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Abstract 



^f) • We improve the McLerran-Venugopalan model for the gluon distribution 

functions in very large nuclei by imposing the condition that the nucleons 
qq , should be color neutral. We find that enforcing color neutrality cures the 

infrared divergences in the transverse coordinates which are present in the 
McLerran-Venugopalan model. Since we obtain well-defined expressions for 
&\ • the distribution functions, we are able to draw unambiguous conclusions about 

various features of the model. In particular, we show that the gluon distri- 
bution functions in the absence of quantum corrections behave as l/x F to 
all orders in the coupling constant. Furthermore, our distribution functions 
exhibit saturation at small transverse momenta. The normalization of the 
distribution function we obtain is not arbitrary but specified in terms of the 
nucleon structure. We derive a sum rule for the integral of the gluon dis- 
tribution function over transverse momenta, and show that the non-Abelian 
contributions serve only to modify the shape of the transverse momentum 
distribution. We obtain a relatively simple expression for the mean value of 
the transverse momentum-squared. The connection between the McLerran- 
Venugopalan model and the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi equa- 
tion is discussed quantitatively. Finally, we illustrate our results in terms of 
a simple nuclear model due to Kovchegov. 

24.85,+p, 12.38.Bx 
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I. INTRODUCTION 



Several years ago, McLerran and Venugopalan JIJ recognized the existence of a regime in 
which the gluon structure functions ought to be calculable within quantum chromodynamics 
(QCD). Their idea is to recognize that in large nuclei at small enough values of the longi- 
tudinal momentum fraction x F , the density of partons per unit area is large. Since a large 
number of charges are contributing, it is expected that classical methods should apply, i.e. 
the vector potential is computed from the classical Yang-Mills equations, and quantum cor- 
relation functions are approximated by ensemble averages. This picture was developed more 
fully in subsequent work |p2|-Q , and is somewhat related to the approach of Mueller p|-p!0| . 

In the region of very small x F , the quantum corrections to the distribution functions 
calculated in the McLerran- Venugopalan model become large @,[TT|]. In particular, these 
corrections are enhanced by powers of a s ln(l/a; F ), implying that the classical approximation 
being made at lowest order is breaking down in this region. This observation has lead to the 
development of a renormalization group procedure [12| whereby these corrections are taken 
into account by incorporating gluons with large values of x F into the charge density. This 
approach, which uses the results of the McLerran- Venugopalan model as its input, has been 



subsequently developed in Refs. |12|-|T8 



The focus of this paper is not on the very small x F region, but rather on dealing 
with the infrared divergence in the transverse coordinates which appears in the McLerran- 
Venugopalan model. The vector potential two-point correlation functions obtained in 
Ref. |12|] are highly infrared divergent, going like (x 2 ) x at large distances. On one hand, the 
region where this poor behavior manifests itself is confined to separations x ^> Aqq D . This 
is firmly in the non-perturbative regime where one does not really believe the results of a 
perturbative calculation anyway. However, since the gluon number density is derived from 
the Fourier transform of this correlation function, the presence of such a strong divergence 
makes it difficult to obtain more than qualitative results from the theory: one must resort 
to cutting off the transverse spatial integrations at Aq^ d |]19 



In this paper we will demonstrate that the infrared divergence described above is closely 
related to the issue of color neutrality. Physically, nucleons appear to be colorless when 
viewed at length scales which greatly exceed Aq^ d . Consequently, we expect that two- 
point correlation functions should die off rapidly beyond distance scales of about Aqq D , 
because of confinement. The quarks in two widely separated nucleons should not feel each 
others presence. Thus, the fundamental requirement that the nucleons be color neutral may 
be rephrased as a restriction on the allowable charge-density correlation functions T>(x). 
Although the McLerran- Venugopalan model employs an ensemble of charge distributions 
which respects the fact that the average charge should vanish, the two-point correlator 
computed from this ensemble is inconsistent with color-neutral nucleons. We claim that the 
color neutrality condition should be viewed as of primary importance, and that we should 
modify the ensemble of charge distributions to satisfy it. By doing so, we incorporate 
one of the major effects brought about by confinement, parameterizing it in terms of the 
precise form chosen for D[x). Enforcing color neutrality cures the infrared divergence in 
the McLerran- Venugopalan model by introducing a new scale into the problem, namely 
the minimum length scale for which color neutrality holds. This outcome is somewhat 
similar to the one considered in Refs. P0,pT[; however, neither of these papers examine the 
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consequences of enforcing color neutrality in great detail. 

Given a model which incorporates the color neutrality of the nucleons at some scale, we 
obtain infrared finite correlation functions. Because our results are well-defined everywhere, 
we are able to draw firm conclusions about various properties of the gluon distribution func- 
tions, in spite of their highly nonlinear form. We are able to demonstrate that the gluon 
distribution functions in the McLerran-Venugopalan model are proportional to l/x F , inde- 
pendent of the distribution of charge in the longitudinal coordinate. To evade this conclusion 
apparently requires charge-density correlation functions which have a non-trivial longitudi- 
nal structure. Although we have modified the form of T>(x) from the one employed in 



Ref. ||12|| , we still obtain correlation functions which are consistent with asymptotic freedom. 
That is, in the ultraviolet limit, the non-Abelian terms drop out and we obtain a gluon 
distribution function which goes like 1/q 2 . At small values of the transverse momentum 
q, our gluon distribution functions either approach a finite constant or diverge mildly (like 
\nq 2 ), depending on the long distance behavior of T>(x). The density of gluons at small 
q grows more slowly with increasing nuclear size than the number of nucleons, consistent 
with saturation [|1(],^-|24| . We also derive a transverse momentum sum rule, which shows 
that the non-Abelian corrections only serve to transfer gluons (at fixed x F ) from one value of 
transverse momentum to another: the total number of gluons is not altered. This is reflected 
in the fact that the structure functions we obtain are simply proportional to the number of 
nucleons: all non-linear dependence on the amount of charge present is suppressed at large 
values of the momentum transfer. As an application of the sum rule, we make contact with 
the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi equation |25|j26|| , demonstrating the plausi- 
bility of the McLerran-Venugopalan model. In addition, we are able to compute a fairly 
simple expression for the mean transverse momentum-squared of the gluons which is accu- 
rate to order unity: that is, the logarithms that arise in our calculation are all evaluated 
at a scale which is completely calculable given a nucleon model. We find that the leading 
Q 2 behavior of this quantity is linear in the number of nucleons, whereas the logarithmic 
corrections are quadratic in the charge per unit transverse area. 

The remainder of this paper is organized as follows. In Sec. f|, we review the McLerran- 
Venugopalan model as described in Refs. ]I]-|,|T^]. We allow for a generalized dependence 
of the charge- density correlator on the transverse coordinates in order to set the stage for 
the subsequent discussion. In Sec. [TT1] we describe the conditions which the charge-density 
correlator must satisfy in order to be consistent with nucleons which are color neutral. We 
examine the properties of the resulting gluon number density in Sec. [TV| If the charge- 
density correlator is Gaussian and local in the longitudinal coordinate, we show that the 
gluon number density is always proportional to l/x F . We derive a sum rule for the integral 
of the gluon number density over the transverse momenta, and apply it to the computation 
of the gluon distribution function at large values of the momentum transfer. Our result is 
closely related to the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi equation. Sec. [TV] concludes 
with a discussion of the mean value of the transverse momentum-squared obtained by our 
approach. We illustrate our results in Sec. [V] within the context of a specific nuclear model by 
Kovchegov [pTfl . Finally, Sec. [Vl| contains our conclusions. A discussion of our notation and 
conventions may be found in Appendix |A]. The details of how to evaluate various integrals 
arising in the text are contained in Appendices ||, y, and |D|. 
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II. REVIEW OF THE MCLERRAN-VENUGOPALAN MODEL 



A. General Considerations 



The goal of the McLerran-Venugopalan (MV) model is to compute the gluon distribu- 
tion function at small values of the longitudinal momentum fraction x F = q^i uon / Qtudeon 
for a very large nucleus. We begin our discussion with the assumptions which form the 
basis of the MV model, and outline the steps in the computation of the gluon distribution 
function. Having described the ingredients which underlie the calculations, we will explain 
the conditions which are required for this picture to be applicable. 

Consider a nucleus moving along the z axis at nearly the speed of light. The central 
assertion of the MV model is that, in a specific kinematic regime, we may determine the 
gluon distribution function of the nucleus as follows. First, we are told to compute the 
classical vector potential generated by the valence quarks. Since the valence quarks appear 
as a Lorentz-contracted pancake-shaped distribution of color charge, the current appearing 
in the classical Yang-Mills equations takes the form (see Appendix |A] for a summary of our 
notation and conventions) 

J + (x) = p(x~ , x), 
J (x) = 0, 

j\x) = 0. (2.1) 



There is no x + dependence in Eq. Q2.1| ) in the limit that the source (i.e. the nucleus) is 
moving at the speed of light. Because the intuitive parton model is realized in the light-cone 
gauge p8|-|3T|, we should write the solution in that gauge, i.e. A + = 0. 



By definition, the gluon number density is just the gluon number operator evaluated in 
the (quantum) state of interest. This may be expressed in terms of the two-point vector 
potential correlation function as |23|j29|| 



dN q + f 00 _ f 2 f°° 



dq + d 2 q 47r 3 J-00 



[°°dx- f d 2 x Hdx'- f d 2 x'e- ig+{x ~ x ' ) e iq < x - x ' ) (A^{x-,x)A^(x'-,x')). 

J— 00 J J —00 J 

(2.2) 



In the MV model, the average over quantum fluctuations implied by the angled brackets 
on the right hand side of (|2.2|) is replaced by a classical average over a suitably weighted 
ensemble of charge densities. The ensemble should consist of all of the physically realizable 
states of the quarks (in color singlet combinations) within the nucleons and of the nucleons 
within the nucleus. Then, what we require is an expression for the vector potential in terms 
of the charge density plus the probability distribution for determining how likely a given 
charge configuation is. In the MV model, ensemble averages are computed according to the 
Gaussian weight 

y^p]exp j- J dx~d 2 x J dx'~ d 2 x' Tr \p(x~ , x)p(x'~ , x'j\ ^ 2 ~^ - V' 1 (x - x') j . (2.3) 

This weight implies that the charge density has an average value of zero. There are, however, 
fluctuations about this average value, the size of which are governed by fi 2 (x~), which is 
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related to the average charge density squared per unit area per unit thickness.[] The precise 
statement of the relationship appears in Sec. |T|. The longitudinal dependence of Eq. ( j2.3|) is 
purely local in x~ , reflecting the fact that the large number of quarks which contribute to the 
charge per unit area at a given transverse position typically come from different nucleons: 
hence, they should be uncorrelated. The transverse space dependence of the MV model as 
described in Refs. ||-g,|rj is 

V(x-x') =5 2 (x-x'). (2.4) 
Ultimately, we intend to modify this dependence. Therefore, we will write our expressions in 



terms of T> and insert (|2.4j ) only when we wish to specifically discuss the results of Ref. |L2 

Thus, given the solution for the vector potential in terms of the charge density and an 
ansatz for the charge density correlation function, we may compute an approximation to the 
gluon distribution function. We now describe the conditions necessary for this procedure to 
yield a good approximation to the result which would be obtained in a full QCD treatment. 

First, we note that the charge distribution for the nucleus is taken to be a Lorentz- 
contracted pancake moving along the z axis at (nearly) the speed of light. We wish to be 
in the regime where the partons do not resolve its longitudinal structure. This is only true 
for partons with very small longitudinal momentum fractions. To determine how small the 
momentum should be, we note that the corresponding length scale x~ ~ V?g/«on should 
be larger than the thickness of the Lorentz-contracted nucleus. This thickness is of order 
y4 1 / 3 a/7, where A is the number of nucleons, a the radius of a single nucleon, and 7 is the 
(common) boost factor of the nucleons which comprise the nucleus. Writing the mass of the 
nucleon to be m, we see then that the condition on the momentum fraction reads 

n + 4-1/3 
x F = -#2=- < . (2.5) 

^■inucleon ma 



Gluons with longitudinal momenta satisfying (|2.5|) cannot resolve distances shorter than the 
thickness of the Lorentz-contracted nucleus they see. 

An immediate consequence of being unable to resolve nucleons with different values of 
x~ is that we may view the nucleons as being piled up onto effectively the same value of 
x~ . Thus, each unit of transverse area will, on average, contain a large number of valence 
quarks. Since there are a large number of quarks contributing, the total contribution will 
be in a large representation of the gauge group, allowing us to treat the source in the the 
classical approximation. The large number of quarks also allows us to justify the use of the 
Gaussian weight ( |2.3|) via the central limit theorem. The number of valence quarks per unit 
area in a flattened nucleus of radius R is 3A/nR 2 ~ 3A 1 / 3 /ira 2 , since the nuclear radius 
scales with the number of nucleons as R ~ A l l 3 a. Gluons with a transverse momentum q 2 



x The authors of Ref. |T^] choose to replace the spatial variable x~ with what they call the space- 
time rapidity y, defined by y = yo + ]h(xq /x~). Thus, they write ^ 2 (y) and refer to it as the 
average charge density squared per unit area per unit rapidity. We will stick to the description in 
terms of x~ in this paper. 
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see a transverse area of about vr/q 2 . Thus, to maintain the condition that there is a large 
number of quarks in the region being probed, we should only look at transverse momenta 
such that 

A|cd S q 2 £ ^a-\ (2.6) 

that is, q should be soft enough to correspond to a large enough area, but not so soft that 
the coupling becomes strong. A second factor limiting the maximum q 2 for which this 
treatment is valid is the observation that if the transverse momentum of the emitted gluons 
is too large, then the eikonal approximation we are making breaks down. That is, the no- 
recoil approximation that allowed us to treat the valence quarks as being localized on the 
light-cone fails. 

The bottom line is that in large enough nuclei for small enough values of the longitudinal 
momentum carried by the parton and moderately large tranverse momenta, the number of 
quarks per unit area becomes large. In this regime, we expect that we may compute the 
gluon distribution function at lowest order from the classical solution to the Yang-Mills 
equations. 



B. Non-Abelian Weisacker- Williams Fields 



As outlined above, the MV model makes use of the classical solution to the Yang-Mills 
equations in the presence of a source of the form ( |2.1|) , with J + (x) = p(x~,x) the only 
nonzero component, independent of the light-cone time x + . Ultimately, we will take the de- 
pendence on the light-cone distance x~ to be a delta function, corresponding to the pancake- 
shaped distribution described above. Initially, however, we will work with a finite spread in 
x~ (a multi-layer pancake), to avoid the ambiguities that would otherwise arise in the com- 
mutator terms of the Yang-Mills equations \12\. Furthermore, certain quantum corrections 



may be included by a renormalization group treatment P"2"|-rrS|j which generates a specific x~ 
dependence, reflecting the presence of the Yukawa cloud of quark pairs and gluons around 
the original source. So, it behooves us to make provision for such a dependence. 

The simplest route to the solution for the vector potential in the light-cone gauge for 
a current of the form ( |2.1| ) is to begin by finding the solution in the covariant gauge, and 
subsequently transform into the light-cone gauge @J2^|. Therefore, we begin in the covariant 
gauge, d^A^i^x) = 0, where the field equations ( |A3| ) become 



d A v - 2i 



A»,d,A» 



+ i 



A„d u A» 



A 



A^,A L 



J" 



(2.7) 



We recall that in electrodynamics, the /zth component of the vector potential is generated 
by the /zth component of the current (in the covariant gauge). Analogously, we attempt to 
form a solution where the only non-zero component is A + . Imposing this condition, we find 
that this is possible, provided that 



d 2 A + - 2i 



A + , d+A~' 



P- 



(2.8) 



Since A and A 1 vanish, the condition that we are in covariant gauge reduces to d + A + 
Hence A + is independent of x + and Eq. (|2.8|) reduces to 



0. 



6 



'V 2 A + (x ,x) = p(x ,x) 



(2.9) 



Notice that Eq. (|2.9|) is local in x . This is a consequence of the x + -independence of the 



vector potential, which reduces d 2 = —2d + d + V 2 to just V 2 . The solution of (|2.9|) is 



A + (x ,x) = J d 2 x' G(x - x')p(x ,x') 



where 



t\2 



(2.10) 



(2.11) 



The arbitrary length scale A which appears in Q2.ll ) reflects the lack of a scale in the classical 
theory. This same lack of scale leads to the infrared divergence in perturbation theory. Note 
that in the covariant gauge, the fact that J and A are color matrices has introduced no real 
complications into the solution: the vector potential is still linear in the sources. 

The non-Abelian features make themselves felt when we gauge transform our solution to 
the light-cone gauge. We write the gauge transformation as 

A"(x) = XJ(x)A^(x)XJ- 1 (x) - zf^U^Jtr 1 ^). (2.12) 
Asking that the new gauge be the light-cone gauge, A + = 0, we find that U(x) should satisfy 



z[a_U(x)] = U(x)A + (x). 
The solution to Eq. Q2.13Q is the path-ordered exponential 



U(x) = "Pexp 



(2.13) 



(2.14) 



where the sense of the path-ordering V is such that 
i=i 



U(rt) = l + f;^ f d£ f^d^ ••• fi-'dfr A + {$-,x)--'A + (£,x)A + (£,x), (2.15) 

_•_ 1 J — oo J —CO J —oo 



i.e. the factors are arranged such that x increases as the factors of A + are read from left 
to right. Using ( f2.12j ) and (2.14) to obtain the light-cone gauge expressions yields 



A~(x) = 



(2.16) 



and 



AHx) 



d j XJ(x) 



IT 1 ! 



x) 



dc u(r, x) Wa + (c, aolu-^r , x). 



(2.17) 



The Abelian limit may be obtained from Eq. ( p,17| ) by allowing all quantities to commute. 
Equivalently, we may take only those contributions which are linear in the source. 
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It is instructive to recall the features of the solution in electrodynamics for a point charge 
moving along the x 3 axis at the speed of light. In this situation, Eq. ( 2.171 ) reduces to 



A'(s) = ^eoo^. (2.18) 

The only non-vanishing components of the field tensor produced by the vector potential 
in ( P1SD are 

F"(x) = -±8(x-)^. (2.19) 

These components of the field tensor correspond to the transverse electric and magnetic 
fields: the fact that = indicates that these fields have equal strength. There are no 
longitudinal fields generated in this limit. An observer sitting at some fixed position (6, x 3 ) 
would see no fields except at the instant when the charge made its closest approach. At this 
time a 5-function pulse of transverse and mutually perpendicular electric and magnetic fields 
would be seen. The magnitude of the pulse would be proportional to 1/6. The (unobservable) 
vector potential, which was zero before the passage of the charge, takes on a non-zero value 
for all times afterward. Since the fields vanish at these times, the late-time vector potential 
may be thought of as some particular gauge transformation of the vacuum. 

The overall features of the above description continue to hold when we switch to QCD, 
although the details are slightly altered by the presence of the path-ordered exponentials, 
which introduce a nonlinear dependence on the source. The net effect of these factors is 
to color-rotate the source in a complicated fashion. Nevertheless, the chromoelectric and 
chromomagnetic fields are non-zero only at the instant of closest approach by the charge, 
and the vector potential switches from one gauge transform of the vacuum to a different one 
at this instant. 



C. Structure of the Correlation Function 



Having obtained the solution for A in terms of p via Eqs. ( |2.10 ), ( 2.14[) and (|2.17 ), we now 
turn to the question of actually performing the color averaging. The Gaussian weight ( p.3|) 
argued for above leads to the two-point charge density correlator 



(p a (x~, x)p b (x'-, x')) = 5 ab p 2 {x-) S(x~ - x'-) V(x - x 1 ). 



^2.20) 



Thus, to compute the correlation function (A%(x , x)AAx' , x')), we expand A in terms of 



p with the help of Eqs. (|2.10|) , (|2.14| ), and ( [2.17|) , and perform all possible contractions on 
pairs of p's using ( |2.2(J| ). At the end of this process, we resum the series. The details of this 
calculation may be found in Ref. [1^]. Here, we shall simply quote the result: 



(At(x-,x)A b j (x , ',x')) 



5 ab did'jLix - x') 
i\L L(x - x') 



exp 



N cX (x-,x'~)L( 



x 



X 



(2.21) 



Eq. (|2.21|) depends on two new functions. The first, x( x i x ' )> is the total charge-squared 
per unit area at positions to the left of the leftmost of (x~,x'~): 



8 



X(x ,x' ) 



min(a; ,x' ) 



(2.22) 



The appearance of such an expression is not surprising when we recall that the value of the 
vector potential depends on whether or not the charge has yet reached its point of closest 
approach. What xi x ~ \ x '~) measures is the amount of charge in those layers of the source 
which have already passed both of the points which we are comparing. Although the range 
of integration extends to x~ = — oo, in practice this is cut off by the form of fi 2 (x~), which 
for a pancake-shaped charge distribution, should be non-zero only in a relatively small range 
around the position of the nucleus along the x~ axis. 

The second new function appearing in ( |2.21| ) is given by 



L(x 



x 



d 2 Z / d 2 $,'V{i - O \G(x - i)G{x' - O - \G{x - i)G{x - 



\G{x' - i)G{x' - £')] . (2.23) 



Note that this function vanishes when x = x' . Thus, at very short distances the nonlinear 
terms in ( |2.21| ) drop out and the behavior of the correlation function is the same as in the 
Abelian case: 



lim (Af{x~, x)A b Ax'-, x')) = 5 ab X {x-,x'~) d.&Mx - x') 

|£C— a;'| — >0 I J 



(2.24) 



In this limit, all of the dependence on the transverse coordinates is contained in the derivative 
of L: 

di%L{x -x')=j d 2 i J d 2 (-'V{i - t')[diG(x - ii)}ldjG(x> - O] 

= 5 ij C(x-x') + d^x-x'). (2.25) 

To uniquely specify the decomposition appearing on the second line of ( ^.25[ ) we take to 
be traceless. Hence, the only piece of ( 2.25 ) which contributes to the gluon number density 
is 5ijC(x — x'). 



D. Result Assuming Purely Local Charge Density Correlations 



As mentioned earlier, the results presented in Ref. |L2] utilize a purely local correlator for 
the charge density correlation function, i.e. Eq. (|2.4| ). Inserting such a ^-function dependence 
into ( |2.23| ) produces 

L(x - x') = j(x - x') - 7(0) (2.26) 

where, formally,! 



2 The authors of Ref. 12] write 7(2;) = ■^5 2 (x), whereas we obtain the convolution of ^ I 8 2 (x) 
with itself. To demonstrate the equivalence of these two quantities requires that we employ a 
consistent regularization scheme for the infrared divergences, for example by doing the entire 
computation in 2 + e dimensions. 
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7 (aj -x') = J d 2 £G(x - £)G{x' - £). 



(2.27) 



This function contains a quadratic infrared divergence, on top of the arbitrary length scale 
introduced by the logarithmic divergence of G{x — x') itself. Fortunately, it is only the 
combination 'j(x-x') —7(0) which appears in the correlation function. Thus, the quadratic 
part of the divergence drops out. There is still a logarithmic dependence on the infrared 
cutoff. The authors of Ref. |L2| takeP] 



7(a; - x') - 7(0) = —\x - x'\ 2 ln( \x - x'\ 2 A 2 



IQty 



QCD 



^2.28) 



where they have chosen on physical grounds to set this cutoff to Aqq D . The derivative 
of ( gggp is 



di&Mx - x) 



_ ^ 

8tt 
1 

+ — 

47T 



ln( \x 



-r'l 2 A 2 
x I iv QCD 



+ 2 



-6in — 



x — x )Ax — X 



\x 



X 



/|2 



(2.29) 



Note that the terms on the second line of ( 2.29| ) are traceless, and do not contribute to the 
gluon number density (|2.2j). 

Because 7(2;— x') is formally the inverse of V 4 , when we form the trace did' i r y{x — x') that 



appears in ( |2.21| ), we should obtain the inverse of V 2 [i.e. Eq. ( |2.11| )]. Eq. ( |2.29| ) satisfies 



this expectation, since the extra "+2" appearing in the first line may be absorbed into the 
arbitrary scale. What this means is that in a strict evaluation of Eq. fl2.21|) , if we want the 



numerator factor did'p(x — x') to consist of a single logarithmic term, then that logarithm 
must contain a different scale from the one used in L(x — x') itself. On the other hand, in 
the qualitative treatment of the infrared employed in Ref. ||12|| , both logarithms are assumed 
to have the same scale (~ Aqcd), producing 



(Al(x-,x)A«(x'-,x')) 



4(N 2 



x — x 



l\2 



X 



X 



|2a2 



(Af c /167r) X (^- 



-)\ X - 



^2.30) 



According to Eq. ( |2.2|) , the gluon number density depends on the Fourier transform of ( 2.30 ). 
Unfortunately, there is a problem with the correlation function appearing in Eq. ( ^.30| ): it 
diverges like 



x 



at large distances. Thus, the required Fourier transform does not exist 
for any value of the momentum q! In Ref. [fL9|J , the x — x' integration is cut off at Aqq D , 
on the grounds that the approximations made in obtaining ( |2.30| ) are not valid at large 
distances. While physically reasonable, such an approach can only be applied qualitatively. 

It is apparent from the discussion leading to Eq. ( |2.24| ) that Eq. (|2.30| ) nevertheless 
correctly reproduces the Abelian limit at short enough distances. Explicitly, 



3 This expression for 7 differs from the one given in Eq. (26) of Ref. [12| by a factor of 2. This is 
in agreement with the requirement V 2 7(a;) = ^5 2 {x). 
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lira _(^(i-a)^V",«0) = ~^7Z- X(x-,x'-) \n(\x -x'\ 2 A 2 QCD ). (2.31) 



The Fourier transform of this function goes like 1/q 2 . 

At larger distances, <C \x — x'\ <C Aqq D , and for large enough nuclei, the second 
term of ( |2.30| ) becomes negligible. Under these circumstances, the correlation function has 
a Fourier transform which behaves as lnq 2 . To see how large the nucleus has to be to have 
such a regime, consider the function /(£) = £ It has a minimum at £ = e _1 , independent 
of f2, and a minimum value of e~ n ^ e . Thus, to legitimately drop the second term of (|2.30|) 
we should have 

N Y 

n = A C 9 > e - (2-32) 

This constraint is difficult to satisfy: in fact, for experimentally accessible nuclei, such as 
uranium, we actually have Q < 1. Just to obtain Q ~ e requires x ~ 45Aq CD . For the 
model described in Sec. |VJ this implies a nucleus with of order 10 5 nucleons. This is a 
somewhat larger nucleus than is required simply to make the color charge density large 
enough for the approximations described at the beginning of this section to be reasonable 
[compare Eqs. (|2.6| ) and (|2.32|) 1. So while the authors of Refs. |]-f|,12,|19| point out that 



realistic nuclei are only marginally large enough, one must go to even larger nuclei to justify 
dropping the second term in ( |2.30| ). 

Our point is, that the range of momenta where the picture described at the beginning 



of this section and in Refs. [|T^.[Ty.^] is applicable is very narrow at best. While it is 
unlikely that anything short of a full non-perturbative solution to QCD can produce reliable 
predictions for the region \x — x'\ > Aqq D , we may reasonably hope that if we construct 
a model which is well-behaved in this region, then we will be able to trust it closer to the 
boundary point \x — x'\ = Aq^ d than if we use Eq. fl2.30| ), with its unphysical large-distance 
behavior. In the next section we will show that the source of the difficulties encountered 
in Eq. ( [2.30| ) is the inconsistency of the charge density correlator ( |2.20| ) with the intuitive 
picture of color neutral nucleons. When the correlator is altered to reflect this expectation, 
we find that the resulting gluon distribution is well-behaved and completely infrared finite. 



III. COLOR NEUTRALITY 

As described in the previous section, the MV model employs the Gaussian weight ( |2.3|) for 
performing ensemble averages. Thus, the ensemble average of the charge denisty vanishes [|TJ : 

(p a (x-,x)} = 0. (3.1) 

This says that at any given point, when averaging over all of the different nuclei within the 
ensemble, we see no net charge. Eq. (|3.1|) is gauge-invariant in the sense that if we perform 
the same gauge transformation on the fields present in all of the nuclei, the ensemble average 
remains unchanged. 

Let us consider the situation inside a single hadronic state, which we will denote by 
\h). This state could consist of a single baryon or meson, or it could be some collection of 
hadrons. For such a state, the charge density at a given point is not necessarily zero: 



11 



(h\p a (x-,x)\h)^0. 



(3.2) 



However, if we sum over all of the space occupied by the hadronic state, we should obtain 
zero: 

(h\ J dx~d 2 x p a (x~, x) \h) = 0. (3.3) 

Eq. (|3.3|) says that any physically observable hadronic state is a color singlet. Since this is 
true for all of the hadronic states in our Hilbert space, we may use the polarization trick to 
generalize ( |3.3| ) to 

(h\ J dx~d 2 x p a (x~, x) \ti) = 0. (3.4) 



Since we have phrased Eq. ( |3.4|) in terms of color singlet physical states, it must hold 
independent of the gauge chosen for p a (x~,x). 
Now consider the bilinear combination 

Z = (h\J dx~d 2 x J dx'~d 2 x' p a (x- 1 x)p b (x'- 1 x') \h). (3.5) 



On one hand, by inserting a complete set of hadronic states and applying Eq. (|3.4|) we 
conclude that Z vanishes. Since the ensemble utilized in the MV model consists of a specific 
collection of hadronic states [i.e. nuclei with some fixed number of nucleons), we further 
conclude that the ensemble average of Z vanishes. On the other hand, we may also compute 
the ensemble average of Z with the help of the correlation function ( |2.20| ): 



(Z) = J dx~d 2 x J dx'~d 2 x' (p a {x-,x)p b {x'~,x')) 

/POO f 
d 2 Y. J dx~ fj?(x-) J d 2 AV(A). (3.6) 

The integral over S is just the (non-zero) transverse area of the nucleus. Since pi 2 (x~) is 
positive (by definition), the only way to obtain a vanishing ensemble average is to require 
that T> satisfy 

'd 2 AD(A) = 0. (3.7) 



This relation is not true for the choice ( |2.4| ) employed in Ref. fL2| . Hence we conclude 
that this choice is incompatible with the expectation of color- neutral nucleons. In fact, this 
violation of color-neutrality contributes to the poor infrared behavior of the MV correlation 



function (gjMj ) 



Although the function £>(A) is itself not gauge-invariant (its transformation properties 
are determined by the condition that the Gaussian weight be gauge- invariant), the color- 
neutrality condition (3"77|) was, nonetheless, derived in a gauge-invariant fashion, with the 



help of Eq. (|3.4| ). The constraint represented by Eq. ( p.7|) reflects the physical observation 
that the nucleus does not carry a net color charge. 

The fact that a color neutral charge distribution posseses an intrinsic scale explains why 
the infrared behavior of the gluon number density is improved by enforcing ( |3.7| ). We may 
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see the roots of this result by considering Eq. ( |2.10| ), the expression for the vector potential 
in the covariant gauge. As pointed out earlier, it depends on an arbitrary length scale A, 
because of the infrared divergence in the inverse of V 2 . If we integrate both sides of ( [2.10| ) 
over x~ and take the expectation value for the hadronic state \h) [so that we may apply ( p.3|) 
to the term involving A] we obtain 



{h\ / dx~ A + (x-,x) \h) = (h\ / dx~ / d 2 x'G(x - x') p(x~, x) \h) 

J— oo J — oo J 

= —{h\ / dx~ / d 2 x' \n(x - x') 2 p(x-,x') \h). (3.8) 

4:71 J-OO J 

Instead of an arbitrary scale, the logarithm is supplied with the scale intrinsic to p, that is, 
the scale over which (|3.3| ) is true. 

Therefore, let us employ a charge-density correlation function T> which satisfies ( | . 7| ) . 
Furthermore, let us make the reasonable assumption that T> is rotationally invariant. Under 
these circumstances, the function Lix — x') is completely free of infrared difficulties! In this 
situation, we find that ( p.23|) may be simplified to 

L( x -x') = — [ d 2 A V(A) \(x -x' - A) 2 ln(a5 - x' - A) 2 - A 2 In A 2 1 . (3.9) 

167T J L J 

(see Appendix 0). In addition, the expression for C becomes 

C{x -x') = — f d 2 A D(A) In(aj - x' - A) 2 (3.10) 

In the next section we will study the gluon number density implied by these two equations. 

At various points in our subsequent discussion it will prove convenient to take the form 
of V(A) to be 

T>(A) = 5 2 (A) — C(A). (3.11) 

In Sec. |LV A we will see that to reproduce the Abelian limit at short distances, T>(A) must 
contain a delta function. Likewise, such a delta function is argued for in Refs. Such 
a term is consistent with nucleons that contain pointlike quarks. In addition, we will see 
that such a structure arises naturally in the model examined in Sec. [V[ To satisfy Eq. ( |3.7|) , 
it is clear that the other function appearing in ( |3.11 ) must obey 

J d 2 AC(A) = l. (3.12) 

We conclude this section by quantifying the connection between p 2 (x~) and the total 
color charge per unit area. We begin by defining to be the integral of p 2 {x~) over all 
x~: 



= I dx~ p 2 (x-). (3.13) 

Then, the trace of Eq. fl3TB) reads 
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J dx~d 2 x J dx'-d 2 x' (p a (x-,x)p a (x'-,x')} = (N 2 - l^oovri? 2 J d 2 AV(A). (3.14) 



Imagine that, instead of a color-neutral distribution of charge satisfying (|3.7|) , we consider 
a single quark. Then, it would indeed be true that T>(A) = 5 2 (A), and the charge-squared 
of that quark should be given by 

g 2 C F = (N 2 - yXoo* R 2 (3.15) 

where C F = (N 2 - 1)/2JV C is the Casimir factor for SU(N). Thus, we see that the color 
charge squared per unit area is (N 2 — l)-^, and that the precise meaning of /x 2 (x") is 
actually the charge squared per unit area per unit thickness divided by N 2 — 1. 

IV. GLUON NUMBER DENSITY 

We begin our discussion of the properties of the gluon number density resulting from a 
nucleus consisting of color- neutral nucleons by inserting ( |2.21| ) into the master formula (|2.2| ) . 



Employing sum and difference variables for the transverse integrations we obtain 

dN '^-'fft/rt^^) 



dq+d 2 q 2tt 3 N c J J L(A) 

/OO (-00 _ f r , ~1 

dx~ dx'-e- iq (x ~ x ) exp N cX (x-,x'~)L(A) -1 . (4.1) 

-oo J — oo I L - 1 J 



The result of the S integration appearing in Eq. ( |4.1| ) should be interpreted as nR 2 , the 
transverse area of the Lorentz-contraced nucleus.^ We would like to perform the transverse 
integrations. All of the longitudinal dependence is contained in the function x(x~,x /_ ). 
Recall that within the Gaussian ansatz (|2.3|) , the function \{ x ~ i x '~) is given by 



/min(x ,x' ) 

rfr/i 2 (D- (4.2) 
-oo 

In order to obtain a more useful expression for x, let us define 

X{x-)= [ X c/r/i 2 (D- (4-3) 

J — oo 

Then, we see that Eq. ( |4.2| ) has the form 

x ( x -,x'-) = X{x~) &(x'~ - x-) + X{x'-) e(x~ - x'-). (4.4) 

For convenience, let us assume that the entire nucleus is localized around some positive value 
of x~ . Then, we may insert Eq. ( |4.4|) into (|4.1| ) and do one of the longitudinal integrals in 
each term, with the aid of the identity 



l See Appendix^ for a more careful discussion of this point. 
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/oo 
dxe i{k+i£)x Q(x-y) 
-oo 



i e i(k+ie)y 

k + is 



y>0. 



(4.5) 



The result is 

dN q+ N- 



dq+d 2 q 2vr 3 N c 



2s 



L(A) 



dx e 



-2ex~ 



exp 



N c X(x-)L(A) 



1 



(4.6) 



To get a finite nonvanishirig result, the remaining integral must produce exactly one inverse 
power of e. To see how this comes about, let us integrate the key factors of (|4.6|) by parts: 



2s 



dx e 



-2ex~ 



exp 



c 

-2ex 



N r ,X(x')L(A)\ - 1 

N c X(x-)L(A) 



cxp 



- 1 



+ 



dx~ iV c/ u 2 (ar)L(A)exp N c X(x~) L(A) 



-2ex' 



(q 



:-2 ' 



(4.7) 



Physically, we know that for x~ — > — oo, A"(x _ ) vanishes identically, whereas for x~ — > oo, 
Af(x~) becomes proportional to the total charge squared per unit area.Q Thus, the surface 
term does not contribute. We may let e — > in the remaining integral without mishap, since 
the remaining longitudinal integration will turn out to be finite. Therefore, we find that 



^— = ^£-^ 7rj R 2 ^ / d 2 Ae iq *£(A) rdx-fi 2 (x-)exp N c X(x-)L(A) 
dq + d A q lit 6 q 



The definition ( |4.3| ) for X implies that /x 2 (x ) = dX/dx 
remaining longitudinal integration to obtain the result 



(4.8) 

Hence, we may perform the 



dN 



1 7TR 2 1 



dq + d?q 



N c 2vr 3 



d 2 Ae iq - A 



£(A) 
L(A) 



exp 



KX^UA) 



- 1 



(4.9) 



Remarkably, the gluon number density goes like l/q + independent of the longitudinal charge 
profile! The only nuclear dependence is on the total charge squared per unit area via X^. 
This conclusion depends upon the fact that the Gaussian from chosen in ( |2.3| ) is local in x~ . 
To generate a dependence other than l/q + apparently requires non-trivial correlations in x~ . 
As pointed out in Ref. [12], a complete renormalization group analysis of the longitudinal 
dependence would generate correlations in x~, although including such correlations was 
beyond the scope of that paper. Since the derivation leading to the form of the correlation 
function given in Eq. ( |2.21| ) depends crucially on having locality in x~ , it is clear that the 
inclusion of such correlations is a non-trivial problem. 



5 If X{x ) grows at asymptotically large values of x , then the total charge squared per unit area 
is unbounded and the integral (|4.6|) diverges. 
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A. Transverse Asymptotic Behavior 



The asymptotic behavior at short or long transverse distances of the gluon number 
density ( |4.9| ) is governed by the behavior indicated by Eqs. ( |3.9| ) and ( |3.10|) for the functions 



L(x) and C(x). Consider first the situation at short distances. It is straightforward to 
demonstrate that for small x 

j, x J x 2 , if T>(A) is smooth at A = 0; (A-\n\ 

{X) ~ \x 2 \nx\ if P(A) contains 5 2 (A). 1 Uj 

Likewise, the short-distance form of C is 

r( \ J const., if T>(A) is smooth at A = 0; ( a \ \\ 

L[X) ~ I lna: 2 , if V{A) contains 5 2 (A). [ } 

In either case, it follows that 

lim L(x) C(x) = 0. (4.12) 

x^O 

This observation will be used in the next subsection where we will derive a sum rule for the 
integral of the gluon number density over q. According to ( 4.11J) , if we want to reproduce 



the Abelian limit at large q with a gluon number density which goes like 1/q 2 , T>(A) must 
contain S 2 (A). Otherwise, the necessary logarithmic divergence in the correlation function 
at short distances will not occur, and the gluon number density will fall off faster than in 
an Abelian theory. 

The situation at long distances is a bit more complicated. For a rotationally symmetric 
correlation function T>, the color neutrality condition ( |3.7[ ) may be rewritten as 

/ ciAA£>(A) = 0. (4.13) 
Jo 

Because this integral is convergent, at large distances we have either T>(A) ~ 1/A 2+P where 
p > (and could be non- integral), or else T>(A) oscillates.^ We will not consider the second 
possibility any further, as it corresponds to onion-like nucleons, with successive layers of 
alternating color charge all the way to infinity. 

For positive p, we find from ( |3.9[ ) and ( |3.1U| ) that 



L(x) ~ < 



P [ P) (4.14) 
in x , it p = 2, 



I 16 

where "sgn T>" stands for the sign of T> at large distances. The long-distance form of L is 



6 If T>(A) falls off faster than a power of A, then the infrared behavior is even better than described 
in this section. The model described in Sect. [V| provides two examples of this type of correlator. 
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Combining these expressions into the integrand appearing in ( |4.9| ) we see that 



£(A) 
L(A) 



exp{iV c A' 00 L(A)}-l 




iV c A' 00 sgnI? 
p 2 (2-p) 2 



A 2- P 



16 



ln a A' 



if P 7^ 2; 
if p = 2. 



(4.16) 



Since A'oo is proportional to the charge- squared per unit area, it is positive. Thus, we see 
that whether or not the gluon number density remains bounded will be determined by the 
sign of T>(x) at large distances in addition to the value of p. 

First consider the situation when sgnX? = — 1. Such a correlation function corresponds 
to the notion of a charge which is screened at large distances. If < p < 2, then the 
exponent in ( 4.16Q goes to — oo and we are left with an integrand which falls like 1/A 2 . This 
corresponds to a gluon number density which grows like lnq 2 at low momenta. For p > 2, 
the exponent goes to zero, and we end up with an integrand which falls like 1/A P . The 
Fourier transform of such a function goes like q p ~ 2 at small momenta. Thus, we conclude 
that the gluon number density is constant at small q in this situation. Finally, if p = 2, 
the leading behavior of (|4.16|) gives an integrand which falls like (A In A) -2 , intermediate 
between the two previous cases. Thus, when sgnP = —1, we are guaranteed that the Fourier 
transform appearing in the gluon number density exists. 

On the other hand, when sgnX? = +1, we encounter difficulties for < p < 2. In this 
situation, the integrand blows up at large distances. In order to have a physically sensible 
gluon number density when sgnV = 1, we must have p > 2, in which case we once again 
have an integrand which falls like 1/A P . 

In Ref. [Q, McLerran and Venugopalan argued that for q — > 0, the gluon number density 
should approach a constant. However, as the authors of Ref. |T^] correctly observe, the 
subtleties associated with using a 5-function for the longitudinal dependence source were 
not correctly taken into account in Ref. 0. Consequently, these authors claim that the 
correct q — > dependence is lnq 2 . In our refined treatment, where the nucleons are forced 
to obey color-neutrality, we see that either behavior is possible, depending upon the form 
of the two-point charge density correlation function at large separations. If T>(A) falls off 
faster than 1/A 4 at large distances, then the gluon number density approaches a constant 
as q — > 0, independent of the sign of T> in this region. On the other hand, if T>(A) falls 
off faster than 1/A 2 but more slowly than 1/A 4 , we obtain a gluon number density which 
grows like lnq 2 for q — > 0. In this case, we must have sgnV = — 1 to obtain a physically 
sensible result. 

Our gluon number density possesses two distinct noteworthy features at small transverse 
momenta. The first, as described in the previous paragraph, is the softening of the q de- 
pendence from 1/q 2 . This softening is present for all values of A'oo, and it even occurs in 
an Abelian theory! As we have already noted, the Abelian expression for the correlation 
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function is simply proportional to £(A). Hence, from Eq. (|4.15|) we see that even the long- 
distance Abelian correlator possesses the power-law behavior required to slow the growth 
of the distribution function to less than 1/q 2 at small transverse momenta. The physical 
interpretation of this result is simple: a color-neutral system of size a is an inefficient gener- 
ator of radiation with wavelengths A ^> a. At the end of the next section we will discuss the 
second noteworthy feature of our results: our distribution functions depend on at low q 
in a manner consistent with the picture of gluon recombination at high densities envisioned 



by Gribov, Levin and Ryskin [22 



B. Transverse Momentum Sum Rule 



Now that we have established the conditions under which the Fourier transform appearing 
in (|4.9Q is well-defined, let us see more explicitly how the non-Abelian corrections alter the 
transverse momentum dependence of the lowest order gluon number density. In the Abelian 
limit, Eq. (|4.9|) becomes 



dN 



dq + d?q 



(N 2 - 1) 



lowest 
order 



Xoo-kR 2 1 
2tt 3 ~q+ 



J d 2 Ae iqA £(A) 



(4.17) 



Let us subtract (|4.17 ) from (|4.9| ) and integrate over all transverse momenta. The only 
q-dependence contained in either expression is e* 9 ' A : hence, the q integration produces 
47r 2 5 2 (A). This allows us to trivially perform the A integration, with the result 



j* q 



dN 



dq + d 2 q 



dN 



all 

orders 



dq + d 2 q 
N 2 



lowest 
order 



i r 2 , a a) 

— - hm - ; . ; <> exp 



N c q+ a^o L(A) 



N^LiA) -l-NcX^HA)}. (4.18) 



Since the position-space function L(A) vanishes at zero separation, we may evaluate the 
limit by expanding the exponential, producing 



dN 



dq + d 2 q 



dN 



a ii dq + d 2 q 



N C (N 2 



l)X 2 R 2 — lim L(A)£(A). 



lowest 
order 



(4.19) 



But, according to Eq. ( |4.12 ), the limit appearing on the right hand side of ( fL19j) vanishes. 
Thus, we are left with the transverse momentum sum rule 



h 2 " 



dN 



dq + d 2 q 



dN 



all 
orders 



dq + d 2 q 



(4.20) 



lowest 
order 



Eq. ( |4.20| ) states that the non-Abelian contributions have no effect on the total number of 
gluons: we could have obtained the same number of gluons by ignoring the non-linear terms 
in the vector potential. What these contributions actually do is to move gluons from one 
value of the transverse momentum to another. Thus, the total energy in the gluon field is 
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affected by the non-Abelian terms. Note that this conclusion is independent of the form of 
V(A). 

If we now assume that T>(A) is of the form ( |3.11| ), with the smooth part of the correlation 



function C(A) being a positive monotonic function, we can show from Eqs. ( [B24|) -( |B26| ) 
that L(A) and C(A) are both monotonically decreasing functions of A. Since L(0) = 0, this 
means that L(A) is negative. Comparing Eqs. (|4.9|) and ( |4.17|) we see that the integrand 



for the all orders gluon number density differs from the lowest-order result by the factor 

exp[iV c A' 00 L(A)] - 1 



(4.21) 



which is less than or equal to unity for L(A) < 0. Since L(A) becomes more negative at large 
A, we see that the non-Abelian corrections actually serve to reduce the degree of correlations 
in the infrared, resulting in a depletion in the number of low-q gluons relative to the Abelian 
result. This depletion can be quite drastic. For example, when p < 2, we have an Abelian 
distribution which diverges like 1 / q 2 ~ p at small q whereas the non-Abelian terms moderate 
this to lnq 2 ! On the other hand, in the ultraviolet the number of gluons is unchanged: the 



factor in ( |4 . 2 1| ) is very nearly 1. Hence, our transverse momentum sum rule tells us that 
there must be an enhancement in the number of gluons at intermediate momenta. Indeed, 
we will see this explicitly in the model presented in Sec. |V[ Furthermore, ( [4.21D implies that 
the effect becomes more pronounced as X^ increases. This is consistent with the idea of 



saturation presented in p2| , which is framed in terms of parton recombination. When X^ 
is large, gluons begin to overlap, and can readily merge via the non-Abelian terms in the 
equations of motion, which become more and more important in determining the correlation 
function ( [2.2 1| ) as the density is increased. If we further imagine that each source factor p(x) 
contributes some characteristic amount of transverse momentum, then at the end of the day 
we find fewer gluons with low q, and a corresponding enhancement at intermediate q. 



C. Gluon Structure Functions and the DGLAP Equation 



The usual gluon structure functions resolved at the scale Q 2 , namely gA(x F ,Q 2 ), may 
be obtained by supplying the trivial factors needed to convert q + into x F and integrating 
Eq. ( |4.9|) over transverse momenta less than or equal to Q: 



9a(x f ,Q 2 ) = / d 2 q 

J\a\<Q 



dN 



q \<Q dx F d 2 q 
N 2 -lnR 2 1 

— K 



27T 3 X p J\q\<Q 



— d 2 q d 2 Ae iq - A 



£(A) 
L(A) 



exp 



(4.22) 



Suppose that we consider measuring the gluon distribution function at large Q 2 . Then, we 
may use the transverse momentum sum rule ( |4.20| ) to replace the all-orders expression on 
the right hand side of (|4.22|) by its Abelian counterpart: 



9a{x 1 



, Q 2 ) = — [ d 2 q [ d 2 A e iq -* C( A) (N 2 - V 

X„ J\q\<Q J 



X^ttR 2 
2tt 3 ' 



(4.23) 
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As explained in the discussion surrounding Eq. ( |3.15|) , the combination (N 2 — ^X^R 2 is 
the total charge squared of the nucleus, which for a nucleus containing A nucleons is simply 
3Ag 2 C F . 

Eq. (|4.23| ) is useful because although we are unable to perform the Fourier transform 
appearing in ( |4.22| ) for the all-orders result with an arbitrary correlation function X?(A), we 
are able to do so at lowest order. Recall that C is the convolution of T> with a logarithm 
[see Eq. ( |3.10| )]. Its Fourier transform is just the product of the transforms of these two 
functions. Therefore, Eq. (|4.23|) simplifies to 



9a(x f ,Q 2 



±1 *< 



V{q) 3Ag 2 C F 
2q 2 2tt 3 



(4.24) 



Finally, if the correlation function is rotationally invariant, we may do the angular integration 
to obtain 



g A (x F ,Q 2 ) = 3AC F — / -^P(g)-. 



Q 2 dq 2 



71 



(4.25) 



Thus, for large enough values of Q, we see that the number of gluons at a given value of x F 
scales with the number of nucleons. Thanks to the sum rule, this statement is true to all 
orders in the coupling constant: any non-Abelian corrections to (|4.25|) vanish in the Q 2 — > 00 
limit. 

The expression in ( |4.25| ) is also closely connected to the Dokshitzer-Gribov-Lipatov- 
Altarelli-Parisi (DGLAP) equation p5| , p6| . Denoting the gluon distribution function by 
g(x F ,Q 2 ), the quark distribution function by qi(x F ,Q 2 ) for flavor i, and the antiquark dis- 
tribution function by qi(x F ,Q 2 ), the DGLAP evolution equation for the gluon distribution 
function of a single nucleon reads 



dg(x F ,Q 2 ) = a s {Q 2 
d\nQ 2 



2tt 



1 # 

F £ 



91 



X 



« S (Q 2 



+P< 



<l!l 



a s (Q 2 



i 



(4.26) 



The functions P gq and P gg appearing in Eq. ( |4.26| ) are the usual Altarelli-Parisi splitting 
functions p6| . At lowest order, we set g(£,,Q 2 ) = on the right hand side of ( |L26| ), in 
accordance with the premise of the MV model that the gluons are generated only by the 
valence quarks. Thus, the only splitting function we require is (to leading order) 



P, 



<iq\X) — C F 



1 + (1 



X 



X 



Inserting ( |4.27[ ) into ( f4.26| ) we obtain 



dg{x F ,Q 2 
dhiQ 2 



a s {Q 2 



2n x. 



1 + (1 



x) 



(4.27) 



(4.28) 



limit 



Since the MV model is only supposed to apply at small x F , we take the x F 
of ( f4.28[ ). In this limit, we end up with the sum of the quark and antiquark distribution 
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functions for all flavors integrated over all values of x F : this is just 3 for a single nucleon. 
Hence 

dg(x F ,Q*) a s (Q 2 ) 1 

dlnQ 2 ~ 3Cf ^^- f - (429) 

We now see that the DGLAP equation at small x F , Eq. ( |4.29|) , is very similar to the MV 
result of Eq. ( |4.25| ). Remarkably, the MV gluon distribution function is almost A times the 
DGLAP result for a single nucleon. Until now, we have been ambiguous about the precise 
scale at which we should evaluate the strong coupling. This comparison tells us that we 
ought to let it run, evaluating it at q 2 and keeping it inside the integral in Eq. ( |4.25| ). By 



doing this, we effectively incorporate certain quantum corrections in our otherwise classical 
treatment. The only other difference between Eqs. ( |4.25| ) and ( |4.29| ) is the appearance of 



T>(q) in the MV expression. Since the color neutrality condition Q3.71 ) also tells us that 
T>(0) = 0, we see that the presence of this factor serves to regulate the integral at small 
values of q 2 . On the other hand, as described in the discussion following Eq. (|4.11| ), if we 
believe that the nucleon should contain point-like quarks, then T>(A) contains 5 2 (A) and 
T>(q) ~ 1 at large momenta. 



D. Average Transverse Momentum-Squared 



Because the gluon number density we have obtained in Eq. (|4.9|) is well-defined, we are 
able to explicitly compute the average transverse momentum-squared (q 2 ) associated with 
this distribution. In performing this calculation, we have assumed that 



V(A) = 5 2 (A) 



C(A). 



(4.30) 



That is, we choose a rotationally invariant correlator which is consistent with nucleons that 
contain point-like quarks. Furthermore, we work to next-to-leading order in the momentum 
expansion, since the first non-trivial nuclear dependence enters in at that level. 

As we show in Appendix O, the gluon number density may be cast into the form 



dN N 2 



ttR 



71- 



Q 2 



dx /i 2 



x 



d x 



HQx) 
Qx 



£(x)e 



N c X{x~)L(x) 



(4.31) 



As we have remarked earlier, all of the non-Abelian effects are contained in the exponential 
factor. We should like some way to organize the expansion of ( |4 . 3 1| ) in powers of 1/Q. In 
the limit Q — ► oo, the Bessel function appearing in this expresion tends to a 5-function. 
This suggests that the x integration will be dominated by the small x region. To connect 
an expansion in x with the desired expansion in 1/Q we rescale the integral in (f4 .31 ) via 
x = y/Q: 



dN N 2 



7T 



dx fi 2 (x 



d 2 y^l 

y 



£(|jox P 



N r X(x~)L 



(4.32) 



In Appendix 0, we argue that for a correlation function of the form ( |4.30| ), the functions 
L(x) and C(x) are simply some power series in x plus lnx 2 times some other power series 
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in x.f] In the case of L(x) both series start at x 2 rather than x°. Furthermore, we observe 
that the intrinsic length scale which must appear in the logarithms should be of order the 
nucleon size a. This same length scale serves to make the expansion parameter dimensionless. 



Looking at Eq. (|4.32|) , we see that by expanding everything except for the Bessel function, it 
becomes apparent that the final result will consist of terms of the form [l/(aQ) m ] ln n (aQ) 2 . 
Because the series for L(x) begins at x 2 , the smallest inverse power of aQ which may multiply 
ln n (aQ) 2 is l/(aQ) 2n ~ 2 . Furthermore, if we keep track of the powers of X, we see that each 
additional power of X comes with an additional factor of 1/Q 2 - Thus, the leading order term 
is the same as in an Abelian theory. Since we conclude that the non-Abelian 

corrections are enhanced in very large nuclei: the subleading terms become important if 
NcX^-Q 2 . 

We now proceed directly to our results, deferring the details of this somewhat lengthy 
calculation to Appendix 0. At next-to-leading order, we find that 



dN 



1 Z< 



i - 



N r X n 



AnQ 2 



where 



In (J^j 2 = 2( lE - In 2) + J d 2 A C( A) In A 2 



(4.33) 



(4.34) 



7# is Euler's constant, and Z 2 = 3Ag 2 Cp is the total charge-squared of the nucleus. As 
advertised, the scale appearing in the logarithm does not enter in as an additional parameter 
which must be determined from the data. Instead, given a particular model for implementing 
the color neutrality condition ( |3.7|) , we obtain a specific value for the momentum scale, 
Eq. ( [4.341) . The integral appearing in Eq. (|4.34j) is closely related to the value of C(x), i.e. 



J d 2 AC(A)lnA 2 



lim 



8nC(x) + lnx' 



(4.35) 



As we have already remarked in the previous subsection, at leading order in Q, dN/dx F 
simply proportional to the number of nucleons: all of the non-Abelian contributions are 
suppressed for Q — ► oo. Although at first glance the leading order contribution to Eq. (|4.33|) 
looks very different from the result presented in the previous subsection, the two expressions 
actually do agree. To see this, insert the Fourier integral representation for T>(q) into 
Eq. (|4.24j) , perform the q integration, and subsequently integrate by parts. Calculation of 
the sub-leading term requires the full treatment given in Appendix 0. 



Finally, the next-to- leading order expression we obtain for (q ) reads 



Q 2 -4nC(0) N c X a 
ln(Q/Q ) 2 8tt 



(4.36) 



7 If C(A) is derived from a rotationally invariant function of the vector A, then the odd powers 
of x will be absent in this expansion. 
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We see that at leading order, Eq. (|4.36|) is independent of the size of the nucleus. A mild 
dependence on the modelling of the nucleons enters in through the appearance of C(0). 
The non-Abelian terms make their first appearance at next-to-leading order. As previously, 
the scale of the logarithms is not arbitrary, but rather determined by the form of C(A), 
according to Eq. ( |4.34| ). 



V. A MODEL WITH COLOR-NEUTRAL NUCLEONS 

To illustrate some of the features of our improved treatment, we turn to a specific model 
introduced by Kovchegov [p7[ . In this model, we imagine a nucleus of radius R, containing 
A nucleons of radius a ~ Aqq D . Each "nucleon" is made up of a quark-antiquark pair. The 
view in the rest frame of the nucleus is as follows: The quark and antiquark for the ith 
nucleon is assumed to be located with equal probability anywhere within a distance a of 
the center of the nucleon. The nucleons are assumed to be located with equal probability 
anywhere within a distance R of the center of the nucleus. In this manner, confinement of 
the valence partons within the individual nucleons is ensured. For the sake of comparison, 
we will also present results where the quarks are distributed within the nucleons using a 
Gaussian weight instead of the uniform weight employed in Ref. |27 . 

Kovchegov begins with the laboratory frame charge density 

N 

p(x~, x) = gJ2(T a )(Tn [5(x~ - xr)S 2 (x - x t ) - 5(x~ - x'r)5\x' - ajj)] (5.1) 
i=i 

for the ultrarelativistic nucleus. In Eq. ( |5.1| ) the quark belonging to the ith nucleon is located 
at (x~ , Xi) whereas the ith antiquark is located at {x'f, x^). The T a 's are generators in the 
SU(iV c ) color space, while the Tf's are similar generators in the color space of each nucleon. 
Kovchegov explicitly constructs the light-cone gauge vector potential generated by this dis- 
tribution, and uses that result to derive an effective two-dimensional charge distribution 
which leads to this vector potential, i.e. he writes 

p(x~,x) = 8(x~)p 2 (x) (5.2) 

and computes p 2 {x). From this explicit form for p 2 (x), one may then average over the 
positions of the quarks and nucleons as described above to obtain the two point correlator 
(p 2 (x) p 2 (x')) . This quantity is connected to the three-dimensional correlator of Eq. (|2.20 ) 
via 

(p a (x-, x)p b (x'-, x')) = 6{x~ - x") (p a 2 (x)p b 2 (x')) (5.3) 

In this way, Kovchegov's model specifies the ^-dependence of the charge density correlator. 
The result obtained in Ref. |27| is of the form 



(p a 2 (x)p b 2 (x')} = £ A6 ab [l(x, x') - l(x, x')] . (5.4) 

In Eq. ( |5.4j ), the function T(x, x') is the contribution that results in the averaging when a 
quark overlaps with a quark (or an antiquark overlaps with an antiquark), whereas T(x, x') 
is the contribution that results when a quark overlaps with an antiquark. Explicitly, 
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d 2 r dz f d 2 $, d( 



., rd^rdz r 



4 s 



5 2 \ 



X 



r-£)5 2 (x' -r-i) 



and 



l(x, x') 



d 2 rdz rd 2 £d( rd 2 £'dC,' 



Ivri? 3 



4 S 

g7TGr 



4 S 

g7r cr 



5 2 (x-r-£) 5 2 (x'-r-g) 



(5.5) 



(5.6) 



In Eqs. ( |5.5| ) and (|5.6|) , (r, z) locate the position of the nucleon with respect to the (arbitrary) 
origin while (^) locate the (anti)quark with respect to the center of the nucleon. Note 
that the 5 functions are two-dimensional, whereas the integrations are t/tree-dimensional. If 
we use a Gaussian instead of a uniform distribution, we should make the replacement 



4 ■? 
3™ 



(2vra 2 ) 3 / 2 



exp 



2a 2 



(5.7) 



for the quark positions. 8 A detailed evaluation of all of the integrals relating to the uni- 
form distribution is presented in Appendix [D|. In the Gaussian case, the integrations are 
straightforward, and we shall simply present the results. 

The function defined by Eq. ( |5.5|) is obviously proportional 5 2 (x — x'), the precise x 



dependence postulated by the authors of Ref . |T2| . It is also obvious that Eq. fl5.6p contains 
no such singular behavior: it is, in fact, a smooth function of the coordinates. For physical 
applications where we will be looking at scales much smaller than the nucleon radius a, 
this smooth function may be neglected in comparison to the 5-function terms, as was done 
in Ref. [27fl. However, as we noted in Sec. fll, taking only the 5-function term produces 



a correlator which is incompatible with the idea of color-neutral nucleons. Thus, when we 
look at scales approaching a, it becomes crucial to include the contribution from Eq. Q5.6D . 
Eq. ( |5.4j ) obviously implies that T>(x — x') has the form introduced in Eq. ( p.ll| ). By sorting 
out the various prefactors associated with (|5.3|) , (|5.4|) , and (P3| ), we may determine the value 
of A'oo from the requirement that the two contributions to ( |3. 1 l[ j each have unit integral. 
We find that 



A^ c %R 2 ' 



(5.8) 



Eq. ( p.8p implies a total charge squared of 2Ag 2 Cp, instead of 3Ag 2 Cp as was used in 
Sec. |IV G| , because this model uses "nucleons" which are qq pairs instead of qqq triplets. The 
function C(x — x') which enforces color neutrality on scales near and beyond the nucleon 
radius is given by 

C{x - x') = — I d 2 £ ( d 2 £ Va^eVa 2 - f 2 5 2 {x -x'-£ + g) 

A7T z a b J\£\<a J\?\<a 



8 There is no a priori reason to assume that the parameters a and R appearing in the uniform 
and Gaussian distributions are the same. When we need to make the distinction clear, we will use 
a subscript U on the uniform parameters and a G on the Gaussian parameters. 
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96(1 -X) 



167T0 2 

for a uniform distribution of quarks and nucleons and 

C(x - x') 



(2 + X 2 )Vl -X 2 - X 2 (4 - X 2 ) tanh" 1 Vl - X 2 



Ana 2 



(5.9) 



(5.10) 



for the Gaussian distribution. In Eqs. (|5.9| ) and ( |5.10| ) we have defined the dimensionless 
distance measure X = \x — x'\/(2a) (for details, see Appendix [D]). These two functions have 
been plotted in Fig. |l|. In the uniform case, C(x — x') vanishes identically for \x — x'\ > 2a, 
whereas in the Gaussian case it has, not surprisingly, a Gaussian tail. Both functions are 
finite at the origin. 

Now that we have specified the functional form of P(A), we may return to Eqs. ( p.9| ) 
and (|3.10| ) and evaluate the inputs to the correlation function, L(x — x') and C(x — x'). In 
the uniform case we obtain 



L(x — x') = —— 

Air 



V100 200 800 64 / 

-(! + 2X 2 - -X« + j-X 8 ) tanh- 1 Vl^X 2 " 
V5 4 64 / 



6(1 -X) 



and 
C(x - x') 



20tt 



1-X 2 
32tt 
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— + ln 

20 1 ' 



X s 



(5.11) 



^X 2 - X 4 ) tanh" 1 Vl -X 2 - (14 + X 2 )Vl - X 2 ] 6(1 - X). 

(5.12) 



The asymptotic forms of these functions are 



L(x - x') 



f a 2 X 2 






ln( 




4:71 






a 2 


In 






\ 20tt 
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X < 1; 
X > 1; 



(5.13) 



and 



£{x - x') 



111 



'X 2 \ 7 
+ ^ 



(f) 



(exactly), 



X< 1; 
X > 1. 



(5.14) 



Because C(x — x') vanishes at large distances, the Fourier transform of the correlation 
function ( |2.21| ) constructed from ([5 -11 ) and (|5.12|) will remain finite at zero (transverse) 
momentum. 
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If instead we employ Gaussian distributions for the quarks we obtain 



L( X - X ') = - 



;i + X 2 ) Ei(-X 2 ) + exp(-X 2 ) - ln(X 2 ) - (1 + lE ) 



(5.15) 



and 



C(x-x') = Ei(-X 2 



57T 



(5.16) 



Eqs. ( |5.15| ) and ( |5.16| ) contain the exponential integral function evaluated at negative values 
of the argument. We have employed the conventions of Ref. [j32"1, i.e. 



ek-z) = - r 

J X 



dt 
— e 
t 



>>0). 



(5.17) 



Also appearing in the first of these two equations is Euler's constant je- Asymptotically, 
we have 



a 2 X 2 

4:71 



L(x — x 



ln(X 2 )+ lE -2\, X<1; 
ln(X 2 )-(l + 7£ )l, X>1; 



(5.18) 



and 



C(x - x') 



-^[ln(X 2 )+ 7£ 

1 exp(-X 2 ) 
8^ X~ 2 ' 



X < 1; 
X > 1. 



(5.19) 



In this case, C(x — x') has an exponential tail at large distances, and we end up with a 
correlation function which falls like exp(— X 2 )/[X 2 ln(X 2 )] in the infrared. Again, the gluon 
number density will be finite at zero transverse momentum. 

We have compiled a series of plots (Figs. 0-|^) to aid in the comparison of the two 
models to each other and to the results of Ref. ||12|| . Each version of the correlation function 
effectively depends on two parameters, which we take to be and au, clg, or Aqcd- 
In preparing these plots, we have adjusted these parameters so that the ultraviolet limit 
of ( |2.21| ) is identical for all three cases. Recall that in the ultraviolet limit, appears 



multiplicatively [see Eq. ( 2.24 )]. Thus, matching in the ultraviolet really means matching 
the values of C at short distances. From ( f4.17p , ( |5.14p and ( p.l9j ), we determine that 



a G = 2exp(^ - fja v 
w 0.464a[/ 



(5.20) 



and 



A 



QCD 



1.7/4, 
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lA4.au 1 . (5.21) 



In Fig. (Q), we have plotted the trace of the correlation function ( p. 21 ) in position space 



versus the dimensionless distance X. In drawing these curves, we have assumed that the 
longitudinal coordinates have been fixed at some value such that a 2 x(x~ ,x'~) = 20. Notice 
that the curve for the MV correlator given in Eq. ( |2.30| ) begins to diverge significantly from 
the other two curves near X = 0.3. At X « 0.347, corresponding to \x— x'\ = Aqq D , the MV 
correlator vanishes. Beyond this point, it rockets off to — oo faster than any exponential. In 
contrast, both the uniform and Gaussian curves are well-behaved. The correlation function 
generated from the uniform quark/nucleon distribution vanishes for X > 1. 

Next, we perform the Fourier transform to form the gluon number density, Eq. ( (4.9|) . 



The results are plotted in Fig. (|3|) for a fixed value of the longitudinal momentum q + . In 
order to define the Fourier transform in the MV case, we have done what was suggested in 
Ref. and simply cut off the A integration at Aq^d- All three curves have generally the 



same over-all shape: a plateau at small values of q, a sharp decrease at intermediate values 
of q, and a tail for large values of q. The differing values for q — > may easily understood 
from the position space functions in Fig. |2|. The MV curve has the smallest number of 
zero momentum gluons because of the abrupt cut-off imposed at A = Aqq D . There are 
more zero momentum gluons in the uniform case, since the corresponding position space 
function extends out to \x — x'\ = 2a. Finally, the exponential tail at large distances in the 
Gaussian case produces even more zero momentum gluons. At large momenta, the gluon 
number density is supposed to go like 1/q 2 . To illustrate this transition, we have plotted 
the gluon number density multiplied by (aq) 2 in Fig. |4]. In this figure, the uniform and 
Gaussian curves are virtually indistinguishable. Beyond (aq) 2 = 10 4 , they remain flat. The 
oscillations visible in the MV curve are a result of the sharp cutoff of the Fourier transform 
integral. As (aq) 2 grows, they become more rapid, but decrease in amplitude. 

Next, we illustrate the effects of the non-Abelian contributions to the correlation function. 
For Fig. [5], we generate the ratio 

(A%(x-, x)AZ(x'-,x')) all ordcrs e ^x£(*-aO _ i 



(A?(x-, x)A?(x'~, x')) iowcst order N c xL(x - x') 

From this expression, we see that the relative importance of the non-Abelian terms depends 
on the magnitude of \: if \ — > 0, then this ratio goes to 1. Fig. || plots this ratio for 
various values of a 2 \ ranging from 2 to 40 using the uniform version of Kovchegov's model.0 
We see that the non-Abelian terms have the effect of suppressing the magnitude of the 
correlation functions at large distances relative to the Abelian result. In momentum space, 
this translates into a depletion of low momentum gluons. 

Finally, we present Fig. |(J This is a plot of the all-orders gluon density ( |4.9| ) divided by 
the lowest-order gluon density Q4.17| ) at fixed q + for various values of X^. We see that for 



small values of q, there are fewer gluons as X^ increases. At very large values of q, there 



9 If we employ Eq. fl5,8| ) and set R ^4 1 / 3 a, we find that a typical magnitude for o?x is approxi- 
mately 2.5 for uranium. 
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is no change. On the other hand, at small values of q, this ratio is less than one, signalling 
that in this region the number of gluons is no longer simply proportional to the number of 
nucleons. The amount of suppression increases with increasing X^, consistent with gluon 



recombination scenario envisioned in Ref. 22 . Because of the transverse momentum sum 



rule ( |4.20| ), the area under each of the curves should be equal. Thus, we see a pile up of 
gluons at intermediate values of q. These features of the gluon number density depend only 
on choosing a correlator consistent with color screening at large distances (sgnD(A) = — 1 
for large A). 

We close this section by evaluating the parameters appearing in the expression for the 
average transverse momentum-squared within this model. According to Eq. ( |4.34| ), the value 
of the scale associated with the momentum logarithms may be constructed from Eqs. (j5.14|) 



and (|5TT9| ): 



1 exp(— 7£ + |) cltj 1 « 1.615arr , uniform distribution; 



1 exp(— : y)cl g ~ 0.749a G , Gaussian distribution. 



If we relate ajj to a G in the manner specified by the ultraviolet matching condition ( |5.20| ) 



these two expressions become equal. On the other hand, the value of the smooth part of the 
correlation function at the origin appearing in Eq. ( |4.36| ) in the two models is 



4ttC(0) 



9 



uniform distribution; 

(5.24) 



2a 2 / 

-5-, Gaussian distribution. 

a G 



These two quantities differ by a factor of ^exp(| — j E ) « 1.03 when Eq. (|5.20| ) is applied. 

In Figs. and [5] we compare the approximate results for dN/dx F and (q 2 ) obtained 
in Sec. [IV D| with "exact" curves generated by numerical integration of the full non-linear 



expression for the gluon number density. Since Qo does not depend on whether we use a 
uniform or Gaussian distribution for the quarks, and C(0) only weakly so, we have employed 
the uniform distribution in these plots, since its numerical integral is more convenient to 
set up. Fig. [7| is a plot of x F dN/dx F divided by the total charge-squared. In order to see 
more easily how the next-to-leading order approximation ([4.33|) is performing, we have also 
plotted the ratio of the approximate to exact result in the region where the approximation 
begins to diverge from the full all-orders value. We see that our approximation gives an 
excellent description of the full result for (aQ) 2 as low as about 10 when a 2 A' oc = 2. As 
could have been anticipated by studying Fig. |6], for a 2 A' 0O = 40 our approximation begins to 
break down a bit sooner, at around (aQ 2 ) = 50 or so. 

Fig. H compares the exact and approximate curves for (q 2 ) versus Q 2 . The cases X^ = 2 
and Xqo = 40 differ by so little that they are virtually indistinguishable if presented on the 
same plot. This is a consequence of the leading behavior (q 2 ) oc Q 2 . For a 2 A' 00 = 40 we have 
shown the leading order approximation as well as the complete next-to-leading order curve 
generated from Eq. ( |4.36j ), since this value of X^ is large enough for the difference between 
the two approximations to be visible. We have also plotted the ratio of the approximate 
to exact results in the "interesting" region. The behavior of the atX^ = 40 result may 
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be understood in terms of Fig. ||: for large values of X^, there is a pile-up of gluons at 
intermediate values of q 2 . This effect is purely non-Abelian in nature. Now, the leading 
order approximation is completely independent of any of the non-Abelian contributions to 
the gluon number density. Thus, it is not surprising that the leading order expression 
falls short of the full result at moderate Q 2 . The first sub leading contribution happens to 
overcorrect by a small amount (~ 5%). Presumably the next term in the expansion will be 
negative in this region. In any event, the next-to-leading order expression is good to within 
a few percent for (aQ) 2 as low as about 20 for both values of a 2 A' 00 . 



VI. CONCLUSIONS 

We have improved the McLerran-Venugopalan model by introducing a constraint on the 
charge- density correlation function to ensure that it is consistent with a nucleus which is, as 
a whole, color neutral. We find that imposing color-neutrality eliminates the divergences in 
the gluon number density present in the original MV model at small values of the transverse 
momentum, provided that the transverse part of the charge- density correlation function, 
T>(A) is rotationally invariant and falls off faster than 1/A 4 at large distances. In this 
situation, the gluon number density approaches a constant value as q — > 0. To obtain 
a gluon number density which goes as \nq 2 at small transverse momenta, we must have 
T>(A) fall off more slowly than 1/A 4 . This is permissible only if we impose the additional 
constraint that T>(A) < at large distances. Then we can use a T>(A) which falls off nearly 
as slowly as 1/A 2 . 

Because we have an expression which is mathematically well-behaved, we are able to 
demonstrate that the gluon distribution function within the MV model is proportional to 
l/x F to all orders in the coupling constant, independent of the functional form of /x 2 (x~). 
This conclusion hinges upon the choice of a purely local dependence on the longitudinal 
distance in the charge- density correlator. The inclusion of quantum corrections |l2"Hl8f is 
expected to change this situation. 

We have derived a transverse sum rule, Eq. (|4.20| ), for the gluon number density. This 
sum rule indicates that the total number of gluons may be computed as if the theory were 
purely Abelian: the only effect of the non-Abelian terms is to shift gluons from one value 
of transverse momentum to another. As a consequence of this sum rule, we have shown 
that the gluon structure function in a nucleus at large Q 2 is simply given by A times the 
result of the DGLAP equation for a single nucleon. We have shown that if we employ a 
charge density correlator which is consistent with charge screening at large distances, then 
we have saturation |22[|: as the density of color charge X^ is increased, the number of low 



momentum gluons grows more slowly than the number of nucleons. 

We have also presented relatively simple expressions for dN/dx F [Eq. (|4.33 )1 and the 



mean transverse momentum-squared [Eq. (f06|)] as a function of Q, accurate to order 1/Q 2 
and Q° respectively within this model. We are able to compute the scales of the logarithms 
in terms of a single model-dependent function C(A) which describes the manner in which 
charge neutrality is approached at scales near and beyond the nucleon size. Rather re- 
markably, the complicated non-linear structure of the full gluon distribution function may 
be understood in terms of the Q 2 expansion. That is, only the Abelian terms enter in 
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at leading order in this expansion, with the non-Abelian contributions making their first 
appearance at next-to-leading order. 
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APPENDIX A: NOTATION AND CONVENTIONS 



We begin by denning the light-cone coordinates x^: 

(*V*0 — I — ,-y> 3 

_|_ Ju _1_ dj , < 

(A1 » 

The transverse coordinates x 1 and x 2 form a two-dimensional vector which we denote simply 
as x, without the usual subscript "T" or "_L" to avoid excessive clutter. Our metric has the 
signature (— , +, +, +). Thus, the light-cone dot product reads q^x^ = —q + x~ — q~~x + + q-x, 
and x± = —x T . We will think of the time coordinate, and % cLS the longitudinal 

distance coordinate (z). 

Occasionally it will be helpful to switch to sum and difference coordinates. We will 
employ the convention 

E=|(a: + a;')> A = x-x'. (A2) 

This asymmetric pair of definitions has the desirable property of a unit Jacobian. 
The classical Yang-Mills equations are 

where we have employed matrix form, i.e. J u = T a J ai/ , etc. The T a are the normalized 
Hermitian generators of SU(iV c ) in the fundamental representation, satisfying 2Tr(T a T fe ) = 
5 ab . The covariant derivative is 



(A4) 

and the field strength is defined in terms of the potential as 



= d»A v - d v A» 



(A5) 



These definitions assume that the coupling constant g has been absorbed into J v [see, for 
example, Eq. ( [5.1D 1. Note that this is not the same convention as employed by MV in 
Refs. where they choose to explicitly extract the factor of g from the current. We 

define the charge density p(x~, x) appearing in Eq. (|2.1|) such that the total charge Z of the 
nucleus is given by 

Z = J dx~ d 2 x p(x~ , x) . (A6) 

Because our intuitive picture of the parton model is most simply realized in the light-cone 
gauge [31], we elect to work primarily in that gauge. The light-cone gauge is defined by the 
condition 

A + {x) = 0. (A7) 

However, solution of the classical Yang-Mills equations for the type of source we consider is 
simplest in the covariant (Lorentz) gauge. For clarity, gauge-dependent quantities will be 
written with a tilde in the covariant gauge: e.g. A. The covariant gauge is specified by 

d^{x) = 0. (A8) 
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APPENDIX B: THE FUNCTION L(x - x') 



1. General Color-Neutral Correlator 

In this appendix we present the calculation of L(x — x') for a general correlator V(A). 
Our starting point is Eq. ( |2.23| ) written in terms of the sum and difference variables defined 
as in 



L(x -x') = J d 2 A J d 2 £ V(A) [ G(x - S - ~A)G{x' - £ + §A) 

— \G{x — £ — ~A)G(x — £ + |A) 

- \G{x' - £ - §A)G(a:' - S + |A; 

All three terms of ( pi|) may be evaluated by considering the integral 



S(x, x', A) 



d 2 £ In 



|cc-£-±A| 



16tt 2 J V A 2 

To deal with the logarithms, we insert the identity 

B\ r°° dz 



In 



\x 



S + |A| 2 



(Bl) 



:b2) 



(B3) 



Because we will ultimately insert our result for S into (Bl|), we may drop any terms which 



do not depend on A, thanks to the charge neutrality condition (|3.7p . Thus, the A-dependent 
bits of (IB 21) do not contribute, and we are left with 



S(x, x' A) 



da 



db 



— — / rf 2 S exp -a(S -x + \A) 2 - 6(S -x - \Af . (B4) 



167r 2 Jo a io b 
The XI integration is now Gaussian, and easily performed, yielding 

/ * n 1 r°° da r°° db 1 f . , A s9 , . , , A . 

S{x,x',A) = —- — -—-eM-a(x-±A) 2 -b(x' + lA) 
lo7r JoaJooa + o L 



+ 



[o(g - | A) + b(x' + | A)] 2 
a + b 



The next step is to insert 



1 



00 ./ a + 6 

— oil 

o r V t 



(B5) 



(B6) 



and rescale a — > ra, 6 — > r6. This allows us to simplify the exponent, producing 

S(x,x f ,A) = -±r- r% f°°- rf6(l-a-b)e W {-Ta(l-a)(x-x'-A) 2 }. (B7) 
1d7t jo t z Jo a Jo b L ) 

The b integration is now trivial, and the a integration becomes trivial after the variable 
change 
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T 



U 

ail — a) 



Thus, we obtain 



S(x, x', A) = -L ^ exp{-u(x -x'- A) 2 }, 
lore Jo u } 



which, when inserted into (|B1"|) yields 

1 



Lix - x') 



1,00 du 
16% Jo u 2 



d 2 A V(A) 



u(x— cc'-A) 2 g-«A 2 



(B8) 



(B9) 



(BIO) 



Although the u integration in ( B10|) appears to be divergent, for a rotationally invariant 
correlation function T>, the integral actually is finite. To see this, note that for small u the 
integrand of (piOp looks like 



du 



u 



2(x - x') ■ A - (x - x 



l\2 



T>(A.) + finite terms. 



(Bll) 



The term containing (x — x 1 ) 2 vanishes when we perform the A integration, because of ( p.7|) . 
If T> is rotationally invariant, then the angular part of the A integration will cause the 
{x — x') ■ A term to vanish as well. 

Because the integral over u is finite, we may employ a trick reminiscent of dimensional 
regulation and write 



L(x — x') 



1 



16n 



d 2 A 



du 



o u 



2-e 



V(A) 



g-u(cc-a:'-A) 2 g -uA 2 



(B12) 



where we understand that the limit e — > is to be taken on the right hand side as soon as it 
is safe to do so. The u integration may be performed by using the analytic continuation of 



du 



o u 



K+l 



r-f3u _ e -aul 



-T(l- n)[a K - p K ] 



(B13) 



(Eq. (3.434.1) of Ref. P2|). After applying ( |B13| ) to ( |B12[ ) and performing some algebra we 
obtain 



Ux - x) 



1 1 T{1 



+ 



16n e 1 — e 
1 



d 2 AV(A) 



2 A ■ (x - x') — (x — x 



f\2 



167T 



d 2 A V{A) (x-x'~ A) 2 ln(x - x' - A) 2 - A 2 In A 5 



(B14) 



The pole term is the manifestation of the behavior described in ( [Blip . For a rotationally 
invariant function D(A), it vanishes. Thus, we are left with just 



Ux - x) 



167T 



d 2 A V{ A) \{x - x' - A) 2 ln{x - x' - A) 2 - A 2 In A 5 



(B15) 



Differentiation of (pi5|) to obtain explicit expressions for the two functions defined in 
Eq. ( |2.25p produces 
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C(x -x') = — f d 2 A D(A) ln(x - x' - A) 2 (B16) 

87T J 



and 



r.c-^^/^w ^^^ . (bit) 

Recall that £jj was denned to be the traceless piece of didjL(x — a:'). It does not contribute 
to the gluon number density. 



2. Rotationally Invariant Correlator 



If T>(A) is rotationally invariant, we may perform the angular integration appearing in 
Eq. ( p3 1 5|) . Letting if> be the angle between x and A we have 



L(x) = -±- fdAAD(A) 

107T JO 



2- 



x / dip (x - 2xA cos ^ + A 2 ) \n(x 2 - 2x A cos if) + A ) - A 2 In A (B18) 
jo 1 J 



The non-trivial angular integrals appearing in this expression are 
J dtp ln(x 2 - 2xA cos if) + A 2 ) = j 

and 



27rlnx 2 if x > A; 
27rlnA 2 ifx<A, 



2t 



d0 cos?/' ln(x 2 — 2xAcos , ?/> + A^ 



-2-kA/x ifx>A; 
-2-nx/A if x < A. 



(B19) 



(B20) 



Eqs. (B19) and (|B20| ) force us to split up the radial integration into two pieces: 
1 



L(x) = - I dA A V(A) (x 2 + A 2 ) In x 2 - A 2 In A 2 + 2A 2 
8 Jo J 

dA A£>(A)[lnA 2 + 2 . 
Similarly, we may perform the angular integration appearing in ( B16|) to obtain 

1 rx 1 roo 

C(x) = — / dA A£>(A)lnx 2 - - / dA AT>(A) In A 2 . 
A Jo 4 Jx 



(B21) 



(B22) 



Let us further assume that T>(A) = S 2 (A) — (7(A), where C(A) is smooth at the origin. 
The color-neutrality requirement implies that 



/roo I 
d 2 AC(A) = l, or / dAAC(A) = — . 
jo 2n 

Making this ansatz allows us to rewrite Eqs. ( |B21| ) and ( |B22| ) as 



(B23) 
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L(x) = ~^x 2 + ^-x 2 lnf 4) + \ f X dA AC (A) \(x 2 + A 2 ) lnf + 2(* 2 - A 2 ) 
8tt 16n \^ Z J 8 Jo V x l ) 



(B24) 



and 



where 



C(x) = -I ln(g) + \ jT dA AC7(A) In(^) , (B25) 

\nE 2 = 2n dA AC (A) In A 2 . (B26) 
Jo 



Because of Eq. ( |B23|) , the integral defining the length scale 5 in Eq. ( B26|) will be finite 



for any reasonable function C(A). Furthermore, the value of S will be closely tied to 
the scale implicit in C(A), namely the nucleon size a ~ Aq^ d . Since we have stipulated 
that C(A) be smooth at A = 0, we conclude that Eqs. ( |B24| ) and ( |B25| ) are telling us 
that L(x) and C(x) may be recast in the form of a power series in x plus lnx 2 times a 
power series in x. Furthermore, the power series for C(x) begins at a; , whereas the first 
non-vanishing contribution to the series for L(x) is x 2 . This observation will be useful in 
organizing the computation of the mean transverse momentum-squared associated with our 
gluon distribution. 



APPENDIX C: COMPUTATION OF THE AVERAGE TRANSVERSE 

MOMENTUM-SQUARED 

In this appendix we will outline the computation of the value of (q 2 ), the average trans- 
verse momentum-squared for the gluon distribution (|4.9| ). By definition 

„ 2 dN 



I d 2 qq 

J\a\<Q 



2\ _ J\i\<Q dx F d 2 q , . 

q ' = — " dN • 

d l - 



'\q\<Q dx F d?q 

We will do our computation through next-to-leading order. Since the numerator of (|C1|) 



grows like Q 2 in the limit Q — > oo, this means we work to order Q°. Likewise, since the 
denominator of ( |(J1| ) grows like InQ 2 for large Q, we determine it through order 1/Q 2 . 
For the purposes of this calculation we will assume that T>(A) is of the form 

Z>(A) = 5 2 {A) - C(A), (C2) 

where the function C(A) is smooth at the origin. As discussed in Sec. [IV A|, T>(A) should 



contain S 2 (A) if the results are to agree with an Abelian theory in the ultraviolet. This is 
consistent with having point-like quarks in the nucleons. Beyond specifying that C(A) is 
rotationally invariant and has unit integral [to satisfy (|3.7| )1, we will not make any additional 
assumptions about this function. 
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1. Calculation of (q 2 )dN/dx F to Order Q c 



We begin with the numerator of (JCJ1|) , employing Eq. Q4.8| ) for the integrand, since it has 
a simpler transverse structure than Eq. ([4.9|). Thus, our starting point is 



d 2 x I d 2 qq 2 e tq - x C{x)e N ^ {x ~ )L ^ ) . (C3) 
l«zl<<9 



nR 2 — / dx fi 2 (x 



dx c 



2tt 



If we were to evaluate the q integral as it stands, we would obtain a mildly complicated 
combination of Bessel and Lommel functions. Alternatively, we may replace q 2 e iqx by 
—~V 2 e tqx . Then, the momentum integration appearing in ( |C3| ) is easily performed with the 



help of Ref. [32 



d l qe iqx 

q\<Q 



2nQ 



2 Jl(Qx) 



Qx 



(C4) 



Observe that if the q integration were to extend to infinity instead of being cut off at Q, 
then the integral appearing in Eq. (|U4j) would have produced a delta function instead of a 
Bessel function. Hence, we conclude that 



lim Q 2 ^^l = 2n5 2 (x). 

Q^oo ^ Q X 



(C5) 



Eq. (|C5| ) will serve as a useful diagnostic tool in determining the order in Q of each contribu- 
tion as it is encountered. That is, if the application of Eq. ( p5| ) to a given expression results 
in a finite result, we will conclude that there are no contributions to that expression which 
grow as Q — > oo. On the other hand, if the result of applying Eq. flC5[ ) is divergent, then 
we will have to do the integral involving the Bessel function exactly for finite Q, ensuring 
that we correctly determine not only the contributions which grow as Q — > oo, but also the 
subleading terms as well. 

After inserting (|C4]) into (|C3| ), we integrate by parts to shift V 2 off of the Bessel function. 
Thus, we arrive at 



dN N 2 - 1 



nR 2 — Q z 



doty p, IX 3j p 



The derivatives appearing in Eq. 



C(x)e 



N c X{x-)L(x) 



dx fi 2 (x 

) yield 

V 2 C(x) 

+N c X(x-)C(x) 



Qx 



C(x)e 



N c X(x-)L(x) 



V 2 L(x] 



+2N c X(x' 



dX(x) 



+N 2 X 2 (x-)C(x) d,L(x 



diL(x 

2 



^N c X(x-)L(x) 



■ (C6) 



(C7) 



The next steps are straightforward but tedious. Basically we have to examine each term bit- 
by-bit and do the x integration either after applying Eq. (|U5|) (for the constant contributions) 
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or by doing the actual integral with the Bessel function intact (for the contributions that 
grow as Q — > oo). Rather than go through all of this mathematics in detail, let us only 
mention the highlights. 

The first term in Eq. (|C7[) , up to the ubiquitous exponential factor, is purely Abelian in 
origin. Differentiation of Eq. ( |3.10| ) tells us that 

V 2 C(x) = -\V{x). (C8) 

Thus, the contribution from this term is especially easy to compute, as there is no A 
integration. The singular part of T>(x) leads to a contribution proportional to Q 2 , whereas 
the smooth part produces a constant piece proportional to C(0). 

The second and third terms of ( p7| ) are purely non-Abelian in origin. At lowest order 
(i.e. with the exponential set to unity), they are quadratic in the charge-squared. At the 
end of the day, the longitudinal dependence of these terms reads 

/OO c OO FX 

dx~fi 2 (x~) X(x~) = l dx~ I /i 2 (x~)/i 2 (£~) 
-oo J — oo J — oo 

= (C9) 

Thus, we see the result will be independent of the functional form of fi 2 (x~), as would have 
been manifest from the beginning had we chosen to begin our calculation with Eq. ( |4.9|) 
instead of Eq. 0. 

Because V 2 L(as) = — 2C(x), the second term of ( |C7| ) contains C 2 (x). As a consequence, 
we encounter the transverse integrals 



r x = fd 2 x^^-\n 2 x 2 (CIO) 
J Qx 

and 



T 2 = I d 2 x^^-\nx 2 \n(x-A) 2 . (Cll) 
J Qx 

The first of these two integrals reduces to 

fiyj- /-OO 

T 1 = — dxMQx) \n 2 x. (C12) 
Q Jo 

This integral may be evaluated by observing that Eq. (6.561.14) of Ref. [^2J yields the 
identity 

p iJl( g lK ^(f)-'EjI±iZ|. (ou.) 

Because e is arbitrary, we may expand both sides of Eq. (|C13| ) and equate corresponding 
powers of e to read off the integral of J\(Qx) times a logarithm to an arbitrary positive 
integer power. In the case of immediate interest we have 
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2tt 

~Q 2 



In 



+ 2 7 i 



(C14) 



We deal with T 2 by first writing 



Qx 



ln(Qx) -InQ ln(cc - A) 2 . 



(C15) 



Next, we rescale x = y/Q in the first term to arrive at a form which possesses a transparent 
Q — > oo limit: 



■My) 



r 2 = 2 1 rf 2 j/ 

-InQ 2 J d 2 x 



In ?/ In ( — — A 
1 Q 



Ji(Qx) 

Qx 



\n(x - A) 2 . 



(C16) 



In the first term of (|C16| ), we may neglect y/Q in the logarithm and do the remaining 
integral with the aid of Eq. (6.772.2) of Ref. ||32| . In the second term we simply apply 
The result reads 



r- 27r 



In 



+ 2 7 i 



In A 2 . 



(C17) 



The third term of Eq. 
transverse integral: 



contains among its contributions only one new non-trivial 



d 2 x 



Ji (Qx) x ■ A 



ln(aj - A) 2 . 



Qx x 2 

We obtain the large-Q limit of ( C18|) by rescaling x = y/Q and writing 

2 2y • A „ / 1 



(C18) 



in, A- | 



In A 



QA 2 



(C19) 



Then, (IC18I) becomes 



rr 1 , a2 f > 2 ^i(y) y - A 

Q l iV — -f- 
2 f „ ^(y) (y-A) 2 



(C20) 



The first term of ( |C20|) posseses a vanishing angular integral, while the second term is 
straightforward to evaluate. Thus, 



2tt 



(C21) 



Finally, we remark that the fourth term of Eq. (|C7|) , which consists of contributions 
which are cubic and higher in the charge density squared vanishes in the large-Q limit. 
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Essentially, this comes about because the small-a; behavior of d{L(x) is Xi[l + In as 2 ]. Thus, 
the factor [diL(x)} 2 contains sufficient powers of x in this region to kill off the diverging 
logarithms. 

After working through all of the necessary algebra, we arrive at the surprisingly simple 
result 



dN 



1 



+ 



x F Air 2 
1 Z 2 N c X a 

X F 47T 2 87T 



Q 2 - 4vrC(0) 



\n(^\+2<y E + J d 2 AC(A) In A 

In + 2 7E + J d 2 A C{ A) In A 



(C22) 



where we have denoted the total charge- squared of the nucleus (N 2 — ^X^R 2 by Z 2 . Each 
occurence of InQ 2 in Eq. ( |C22 ) comes with the same combination of 7# and integral over 
C(A), leading to the definition of Qq presented in Eq. Q4.34 ) . 



2. Calculation of dN/dx F to Order 1/Q 2 



We now turn to the denominator of ( pl|) . Previously, we presented Eq. ( 4.24 ), which 
writes this quantity in terms of the momentum space representation of the charge density 
correlation function T>(q). However, since the numerator has been expressed in terms of the 
position space function C(A), it will be useful to do the same for the denominator as well. 
Besides, Eq. ( |4.24|) is correct only to leading order, and we require the 1/Q 2 terms as well. 



To obtain such an expression, we again employ Eq. ( J4.8| ) for the integrand: 

dN N 2 - 1 _„ 1 roo 



ttR 2 — r dx~n 2 {x-) ( d 2 x [ d 2 qe iq - x C{x)e NcX{x ~ )L{x) . (C23) 
dx F Ztt x f J— 00 J 1/ 1 g 1 <<5 

The momentum integration appearing in (|C23| ) may be done by applying Eq. (|C4|) , producing 



™ = ^^.LqI r dx - Axn U x hm ax) e«^-™. (C24) 
dx F n z x F J-00 J Qx 

Since we are working to order 1/Q 2 , it is not enough to simply apply Eq. QC5|) to replace the 



Bessel function by a delta function when it is safe to do so, since this would make order 1/Q 2 
errors. Instead we will have to expand the exponential and apply an assortment of scaling 
arguments to determine the order of each contribution. We will denote the contribution to 
dN/ dx F which contains m powers of N c X(x-)L(x) by dN^/dx F . 
At zeroth order in this expansion we have 

dNi0) _ z2 l ^\f,2 MQx) , 2 f j2A ^f^ MQx) 



dX F 87T 3 



±Q 2 \ [d 2 x^\nx 2 - fd 2 AC(A) [d 2 x^ \n(x - A) 2 ). 
x F [J Qx J J Qx J 

(C25) 
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where we have used Eqs. (|3.10|) and ( |C2| ) to write out C(x). The ^-integration in the first 
term of ( |C25|) may be performed exactly with the aid of Eq. (6.772.2) of Ref. [32|: 



Qx 



. , 27T 



111 



+ 2 7£ 



The second term, which will combine with the logarithm in 
dimensionless, requires a bit more work. To this end we define 



T 4 = Q 2 J d 2 AC{A) J d 2 x 
The rescaling x = wA produces 



MQx) 

Qx 



ln(x - Af 



(C26) 

to make its argument 
(C27) 



% = Q J d 2 A C(A) A In A 2 J d 2 w 



2 Ji{wQA) 



w 



+ Q / d 2 A C(A) A / d 2 w 



Ji(wQA) 



w 



ln(w - Af 



(C28) 



where A is a unit vector in the A direction. We may immediately perform the w integration 
exactly in the first term of (|C28| ). In the second term, we rescale, A = A'/ (wQ): 



T 4 = J d 2 AC{A) In A 2 



+ 



Q- 



d 2 w d 2 A'C 



A 

wQ 



A'MA^—Hw-A) 2 . 



(C29) 



The only angular dependence in the second term is contained in the logarithm. Thus, 
Eq. ( pi9|) may be applied to yield 



T 4 = J d 2 AC{A) In A 2 

27T r°° dlU , n / 



+ ^-\nw 2 d 2 A'C 



Q 2 Ji w 



A' 

wQ 



A'MA). 



(C30) 



Thanks to the theta-function from Eq. (|B19| ), the radial part of the w integration begins at 
w = 1, not w — 0. If we now expand the correlation function in powers of A'/(wQ), the 
w integration will always converge. Likewise, we may always do the A' integration via the 
analytic continuation of Eq. (6.561.14) of Ref. [p^| . Thus, we obtain a well-defined expansion 
of the second term of ( |C30| ) in powers of 1/Q. In particular, the 1/Q 2 term vanishes, since 



rod 

/ dAA 2 Ji(A) = 0. 
Jo 



(C31) 



Thus, only the first term of Eq. ( |C30| ) contributes through order 1/Q 2 . Combining this 
result with Eqs. ( |C25|) and ( |C26|) we arrive at 



dN^ 



x F 4tt 2 



\n(®p\+2 lE + Jd 2 AC(A) In A 5 



(C32) 
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The next term in the expansion of the exponential of Eq. (|C24j) will produce contributions 
of order l/Q 2 . Our overview of this part of the calculation begins with the result of doing 
the expansion, which reads 



dx„ 



Z 2 1 
— NX - /Xi 
q 



2tt 



Q 2 / d 2 x 



MQx) 

Qx 



C(x)L(x). 



(C33) 



Eqs. ( p.9|) and (|3.10| ) when combined with Eq. ( p2|) tell us that there are three kinds of 
contributions to Eq. ( |U33| ). The first is generated by multiplying the 5 function containing 
terms of C(x) and L(x). It may be evaluated exactly with the help of 



/ dx x 2 In 2 x J\ (Qx) = — 
Jo " Q 6 



In 



+ 2 7E 



(C34) 



which was obtained via the same methods used to integrate Eq. (|C12| ). 

The second type of contribution to Eq. ( |C33| ) consists of terms produced by combining 
a 5 function from one factor [C(x) or L(x)] with the smooth part of the other. Generally, 
these contributions are handled in the same fashion as 7^ [Eqs. (|CJ27| ) ( JO30D ] . In particular, 
we extract the l/Q 2 terms via the rescaling x = wA followed by A = A'/(wQ). The only 
complication that arises is the appearance of the integral 



Q_ 

4tt 



d 2 AC(A) A 3 / d 2 w 



w ■ A 



w 



Jt(wQA) ln(wA) 2 InO - A) 



The angular part of the w integration is governed by Eq. ( |B20| ), leading to 

Q 



r 6 = | / gPac(a) a 3 



dwJ x (wQA) \n(wAy 



e(w 



w 2 ea 



w 



(C35) 



(C36) 



We may apply our usual rescaling of A to the first term of (|C36|) to deduce that it is of 
order (lnQ 2 )/Q 4: , and so may be neglected. However, in the second term, we cannot rescale 
A in this manner since the w integration would diverge. Instead, we write w = v/(QA), 
yielding 



1 r rQ A 
% = —j d 2 AC(A) j o dvv 2 (\nv 2 -\nQ 2 )J 1 (v) 



(C37) 



If we let the upper limit of the v integration go to infinity, we see that the integral remains 
finite. Therefore, we conclude that 



Te Q 2 ' 



(C38) 



correct through order l/Q 2 . 

The final type of contribution to 



which we encounter is generated by multiplying 



the smooth parts of both factors, leading to the formidible-looking integral 



■2 MQx) 



d l x 



Qx 



d 2 AC{A) 



x 



A) 2 \n{x- A 2 ) -A 2 In A 2 



d 2 A'C(A') \n(x - A') 2 . 

(C39) 
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By employing the rescalings A = Ax and A' = A'x followed by x = y/Q, it is straightfor- 
ward to show that the leading contribution to ( p39| ) is of order (In 2 Q 2 )/Q 6 - 
At the end of this lengthy procedure we arrive at 



dN^ 1 Z 2 NcXoo 
dx F x F 47r 2 4nQ 2 



\n(^)+2j E + Jd 2 AC(A) In A 2 



(C40) 



When combined with Eq. (|C32|) , we obtain the result (f4.33|) presented in Sec. [IV D 



APPENDIX D: THE INTEGRALS APPEARING IN KOVCHEGOV'S MODEL 
1. The Charge Density Correlation Function 

In Sec. [V|, we wrote the charge density correlation function derived from Kovchegov's 
model []27| in terms of the two integrals X and X, given by Eqs. ( |5.5| ) and ( |5.6| ) respectively. 
In this Appendix, we will explicitly evaluate these expressions assuming that the nuclear 
radius R is much larger than the nucleon radius a. 

The easier of the two integrals to evaluate is X. Because the nucleus and nucleons are 
spherical and we assume a sharp cut-off in the allowed positions of the nucleons and the 
quarks, we have 

I(x,x') = - 62 ^~ X P I d 2 i[ d 2 r VR 2 - rVa 2 - e 5 2 (x - r - £)■ (Dl) 

4 TT 2 R d a S J\£\<a J\r\<R 

Now an exact evaluation of ( pi|) would require that we satisfy the ^-function by setting 
r = x — £ and integrating over only those values of £ which are compatible with this. 
However, since we take R 3> a, we may approximate by setting r 2 w x 2 in the integrand 
and performing the £ integration over the entire disk of radius a. By doing this, we are 
allowing nucleons which are located near the edge of the nucleus to "pop out" by a (small) 
amount of order a. In this approximation, a straightforward evaluation yields 

J(S ' A) = 2^ f-Jt Q ( R2 - J2(A) - (D2) 

where we have switched to the sum and difference variables ( |A2|) . Now Eq. (|D2|) depends 
upon I] whereas the expressions employed in Sec. do not. The S-dependence of Eq. (|D2|) is 
not really surprising, since the pancake obtained by Lorentz-contracting a spherical nucleus 
should be thinner near the edges than near the center. What we call in Eq. (|5.8|) is 
obtained from ( |L)2| ) by first integrating over S: i.e. instead of simply writing nR 2 for the S 
integral, as was done in the discussion following Eq. ( |4.1| ), we perform that integral including 
the geometric factors indicated in (|D2j) : 

|rf 2 SX(S,A) =5 2 (A). (D3) 

Comparing this result with Eqs. (|5.3|) and ( |5.4|) , we conclude that 
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vri? 2 ^ = L?L. (D4) 

The other integral is somewhat more difficult to evaluate. Doing the longitudinal inte- 
grations and making the R^> a approximation described above yields 



( d 2 £ [ Va 2 - £Va 2 - £' 2 6 2 (A - £ + £') 

1-^Q(R 2 -V 2 )C(A) (D5) 



2ixR 2 V R 2 

where C is the function introduced in Eq. ( |5.9|) . Since we have defined C(A) by extracting 
the same S-dependent prefactor which appeared in ( |L)2| ) , at the end of the day we find that 
D(A) = 5 2 (A) - C(A), as written in Eq. (|3~TT| ). 

A direct evaluation of the integral defining C is difficult because of the 5-function con- 
straint: we are instructed to consider the area of overlap between two disks of radius a 
separated by a distance A. This area must be weighted by the product of the distances from 
the two centers to the integration point (see Fig. |9[). Although this sounds simple, the only 
immediate conclusion which may be drawn directly from the integral is that when | A| >2a 
the function vanishes: the disks do not overlap. 

To deal with the situation when | A| < 2a, we begin by computing the Fourier transform 



of C(A): 

C(q) = J d 2 A e^ A C(A) 



9 



, , / d 2 $,[ d 2 i' [ d 2 A e^ A Va 2 - £Va 2 - f 2 5 2 {A -£ + £') (D6) 

4lT z a b J\£\<a J\£'\<a J 

Since the A integration is over all space, the ^-function constraint can be satisfied for all 
values of £ and £': no messy B-functions are introduced. In fact, the entire expression 
factorizes: 

e ^={^Lf^ ,,<v ^} 2 - < D7) 

Performing the angular integration and rescaling the radial integral to unit range we obtain 



2 



C(q) = <3 / dvvVl - v 2 J {aqv) \ . (D8) 



o 



The radial integral may be performed with the help of Ref. |32| (c.f. Eq. (6.567.1)): 

J 3 /2(ag) y 
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Happily, the half-integer Bessel functions are expressible as trigonometric polynomials. 
Thus, we arrive finally at 



C(q) 



(aqY 



sin(ag) — (aq) cos(ag) 



In spite of the apparent high-order singularity at q = 0, C(0) is nevertheless finite: 



lim C(q) 



1. 



(D10) 



(Dll) 



It is amusing to note that the function implied by ( pil|) for the Fourier transform of T>, 
namely 



V(q) = 1 



(aq) 



9 r 

sin(ag) — (aq) cos(ag) 



(D12) 



is exactly the same function that arises in position space for the parallel spin pair correlation 
function in a noninteracting gas of spin-| fermions |J3]. In the Fermi gas, the role of the 
sphere of radius a containing a uniform distribution of quarks is played by the momentum- 
space filling of energy levels up to the Fermi surface. The condition that the nucleons be 
color neutral on large scales leads to the vanishing of T> at zero momentum. The Pauli 
exclusion principle leads to the vanishing of the parallel spin pair correlation function at 
zero separation. 10 

Returning to the issue at hand, to obtain (7(A), we simply invert the Fourier transform. 
The angular integration is straightforward, and leaves us with 



C(A) 



9 



Ana 2 



dqJo{^) 



1 cos 2q 2 sin 2q cos 2q 



+ -r + 



q° q° 



Q 



Q 



(D13) 



Each of the five terms on the right-hand-side of 



is individually divergent. In order to 



perform the q integration term-by-term, we must insert a regulator. Since the final combi- 
nation of terms in ( |D1$ ) is guaranteed to be finite by ( pil ), the result will be insensitive 
to the details of the regulator employed. A convenient means of regulation is to make the 
replacement 



(D14) 



in ( pT3|) . We may then apply Eqs. (6.561.14), (6.699.1), and (6.699.2) of Ref. @ to obtain 



C(A) 



9 



X 2 1 1 



Ana 2 



2 e 1 



+ 



X 4 1 



(1 -e 2 )(l -e 2 /4) 



10 The opposite spin pair correlation function for the Fermi gas is unity (the exclusion principle 
does not affect particles of opposite spin). This CclSG IS cLIlctlcl gous to the non-charge conserving 
correlator of Eq. (2.4). 
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/X 2 y cosvre 

\T" J e(l - 2e)(2 - 2e)(3 - 2e)(4 - 2e) 



x (3 - 2e)(4 - 2e) 2 iq(-l + e, -\ + e; 1 + 2e; X 2 
-4(4 - 2e) 2 F 1 {-1 + e, -| + e; 1 + 2e; X 2 ) 



+4 2 F 1 (-2 + e,-f + e;l + 2e;X 2 



(D15) 



which is valid for A < 2a. In this expression we have defined the dimensionless distance 
X = A/ (2a). The hypergeometric functions may be expanded, and all of the poles in e 
cancelled. The resulting series may then be converted back into a different (generalized) 
hypergeometric function. The result of this tedious algebra is 



C(A) 



9 



64vra 2 



/X 2 \ 

2X 2 (4 - X 2 ) lnf — J + 8 + X 



2 



3 F 2 (1,1,§;3,4;X 2 ) 



(D16) 



Finally, we notice that all but one of the parameters in the hypergeometric function 
appearing in ( pi6[ ) are integers. Hence, we may eliminate this hypergeometric function by 



combining Eqs. (7.512.12) and (9.111) of Ref. [32] to obtain 

„ i , \ r(p)r(g) 

3 F 2 (a, ,c;p,g;,) = r(6)r(f _ 6 y (c)r(g _ c) 

x f dt I' du t b -\l -ty- h ~ l u c - l {l -u) q - c -\l -tuz)- 
Jo Jo 

The required integrations are all straightforward, yielding 



3 F 2 (1,1,§;3,4; z ) = -t(l-t 



(D17) 



z 



z 



lnf -) + 2tanh~ 1 VT^ z 



2 16 8 , , _ > r- 

+- + — -— I 1 + - Vf- : 

z, 



z z 3 z 2 



(D18) 



Inserting ( p!8|) into (|D 1 6|) at last gives us 
9 



^ ^ 16na 2 
which is true for X < 1. 



(2 + X 2 )Vl-X 2 - X 2 (4 - X 2 ) tanh- 1 y/l - X 2 



(D19) 



2. Computation of C(x — x') 

We now turn to the evaluation of C(x — x'). This function may be determined by 
the same procedure used to determine C(A). Since Eq. ( |3.10|) expresses C(x — x') as the 
convolution of T> with a logarithm, we may use ( |3.11| ) and ( |D10| ) to immediately write 
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(aq) 



9 r n 2 
sin(ag) — (aq) cos(ag) 



(D20) 



Once again we encounter an expression with a fake pole at the origin: the q — > limit is 
actually 



lim £(q) = — . 



(D21) 



So, to invert the Fourier transform, we may employ the same trick as was used on Eq. ( P13|) : 
after the angular integration we replace the Bessel function of order zero by one of order 2e 
and integrate term-by-term. After verifying that all of the poles in e do indeed cancel, we 
are left with 



C(x - x>) = ±- [-(1 - X 2 )(8 + 8X 2 - X 4 ) In(^) 



--(1 - X 2 )(169 - 47X 2 + 7X 4 ) 
6 

1 9X 8 

6 64 



(D22) 



This expression is valid only for X < 1. Eq. flD22D displays the precise form predicted for 
C(x — x') at the end of Appendix [B| it is a power series in X 2 plus lnX 2 times another 
series in X 2 . No odd powers of X appear in this series since we employed a rotationally 
invariant correlation function depending upon the vector x — x', as opposed to an arbitrary 
function of \x — x'\. 

The hypergeometric function may be eliminated from ( |D22| ) by combining ( pi 7] ) with 
Eq. (7.512.12) of Ref. |2[: 



4 F 3 (a,6, c,d;p,q,r;z) 



T{p)T(q)T(r) 



T(b)T(p - b)T(c)T(q - c)T(d)T(r - d) 



ds I dt 



du s b -\l- s) p - b -H c -\l-t) 



q—c—l 



xu d - l (l 



u) 



r-d-1. 



stuz)' 



(D23) 



from which we learn that 



4 F 3 (1,1,|,4;3,5,5^) 



64 
Vz 
128 
9^2 



9 8 



z z a 

, 13 14 
1 + r 



In - +2tanh _1 v / T 



224 64 256 1792 

^7 + ^ 



27^ z 2 



9z 4 



(D24) 



The result of inserting ( p24| ) into ( p22| ) is precisely ( 5.12 ) except for the step-function. For 
X > 1, it is easy to combine the integral representation for £ given in Eq. ( |3.10| ) with the 
integral form of C from Eq. (|5.9|) to show that that C(x — x') vanishes in this region. 
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3. Computation of L(x — x') 



Finally, we come to L(x — x'). We should like to begin with the Fourier transform 
function L(q). Note that ( |2.25| ) implies 



f~L(q) 



2C{q). 



(D25) 



Unfortunately, this relation is not sufficient to determine L(q) uniquely: we may add an 
arbitrary multiple of S 2 (q 2 ) to L(q) without contradicting Eq. (|D25| ). According to the 
integral representation of L(x — x'), L(0) vanishes. Therefore, we could adjust this arbitrary 
term so that this is true when we invert the Fourier transform. Operationally, we may simply 
invert the expression implied by ( |D25| j using the same method we have used on the other 
integrals in this Appendix and drop the pole and constant terms. Although this procedure 
lacks the satisfying feeling that comes when one sees the poles cancel automatically, it does 
produce the correct result, as has been verified by a direct numerical evaluation of the 
integral in ( |3.9|) . The bottom line is that for X < 1 we find 



L(x 



x 



a 

4n 



whereas for X > 1 



X 6 X k 



+ 



x z - — + 
3A 2 



In 



128, 
9A 4 5A 6 



+ 



2 

9X 10 
12800 



+ 



(t) 

5A 8 



24 



384 



4 F 3 (1,1,|,4;3,6,6;A 2 ) 



(D26) 



L(x - x') 



20tt 



93 , / 




hln 


m 


20 ^ 





(D27) 



This latter result may be obtained inserting the integral representation ( J5.9Q of C into the 
integral representation Q3.9|) of L. Using the identity ( |D23|) to evaluate the hypergeometric 
function appearing in ( |D26| ) gives: 



4F 3 (l,l,|,4;3,6,6;z) 



20 / 80 640 m 
— 15 + — + — 



9z V 



lnf - ) +2tanh _1 VT 



9z 2 
500 



^ r 982 2984 1488\ , 

25 + —)Vl-z 

z Z A Z A J 

8000 1600 6400 3968 



27z 27z 2 z 3 3z A 3z 5 
Eqs. (|D26|) - (|D28|) combine to reproduce the expression presented in Eq. (|5.11|) 



(D28) 
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FIGURES 
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FIG. 1. The smooth part of the two-point charge density correlation function in Kovchegov's 
model |27j, as given in Eqs. ( |5,9|) (uniform quark/nucleon distribution) and fl5.10|) (Gaussian 
quark/nucleon distribution). The nucleon size parameters a have been chosen according to 
Eq. ( |5.2C| ) so that resulting gluon number densities match in the ultraviolet limit. In the uni- 
form case, C(A) vanishes for A > 2a. In the Gaussian case, there is a Gaussian tail in this 
region. 
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FIG. 2. The trace of the two-point vector potential correlation function in position space. The 
Fourier transform of this function is proportional to the gluon number density [c.f. Eq. ( |2.2j )1. The 
nucleon size parameters a and Aqcd have been chosen according to Eqs. (5.20) and ( |5.21| ) so that 
these functions match in the ultraviolet (\x — x'\ — ► 0) limit. The longitudinal coordinates (x~ , x'~) 
have been fixed at a place where a?x = 20. Plotted are the results of Ref. @ (labelled "MV") 
as well as results using Kovchegov's model [^7) with a uniform or Gaussian distribution of quarks 
and nucleons. 
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FIG. 3. Plot of the gluon number density ( |4,9[ ) at fixed q + . The nucleon size parameters a 
and Aqcd have been chosen according to Eqs. ( 5,20 ) and ( |5.2I ) so that these functions match in 
the ultraviolet limit. The total charge squared per unit area has been set to a 2 Xoc = 20. Plotted 



are the results of Ref. |l2j (labelled "MV") as well as results using Kovchegov's model [27] with a 
uniform or Gaussian distribution of quarks and nucleons. 
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FIG. 4. Plot of (aq) 2 times the gluon number density ( [4.S| ) at fixed q + . This combination has 
been chosen to enhance the high momentum part of the distribution and show the approach to the 
predicted 1/q 2 ultraviolet behavior. The nucleon size parameters a and Aqcd have been chosen 
according to Eqs. ( |5.20| ) and (5.21) so that these functions match in this limit. The total charge 
squared per unit area has been set to a 2 ^^ = 20. Plotted are the results of Ref. El (labelled 



"MV") as well as results using Kovchegov's model [27 with a uniform or Gaussian distribution of 
quarks and nucleons. 
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FIG. 5. The effect of the non-Abelian terms on the trace of the correlation function. Plotted 
is the ratio of the all-orders result to the lowest-order (purely Abelian) result using the uniform 
version of Kovchegov's model [27|. The different curves are for fixed values of the longitudinal 
coordinates (x~,x'~) such that a 2 \ = 2, 5, 10, 20, and 40. 
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FIG. 6. The effect of the non-Abelian terms in momentum space. Plotted is the ratio of the 
all-orders gluon number density to the lowest-order (purely Abelian) gluon number density at fixed 
q + using the uniform version of Kovchegov's model |27|] . The different curves are for a 2 Xoo = 2, 5, 
10, 20, and 40. 
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FIG. 7. Top: Gluon structure function x F g(x p ,Q 2 ) = x F dN/dx F divided by the total color 
charge-squared for a large nucleus. Plotted is a comparison of the next-to- leading order approx- 
imation ( [4.33 ) with the exact result obtained by numerical integration of ( [4.22| ) for two different 
values of Xoo . Bottom: The ratio of the next-to-leading order approximation to the exact result in 
the region where the approximation begins to break down. 
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FIG. 8. Top: Mean transverse momentum-squared as a function of the parton resolution scale 
Q 2 . Plotted is a comparison of the next-to-leading order approximation (|4.36| ) with the exact 
result obtained by numerical integration for a^X^ = 2, and a comparison of the leading order, 
next-to-leading order, and exact results for a 2 Xoo = 40. Bottom: The ratio of the approximate 
results to the exact results in the region where the approximation begins to break down. 
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FIG. 9. The integration region encountered in Eq. ( p5| ) for the smooth part of the correlation 
function in the uniform version of Kovchegov's model p7[ . 
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